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OVERVIEW 


1  . 1  Program  Background 

Linear  acceleration  forces  are  encountered  frequently  in  high-performance 
aircraft  during  rapid  turns  and  pullouts  from  dives-  It  has  long  been  recognized 
that  these  forces  subject  the  pilot  to  a  psycho-physiological  stress  that  could 

■U 

degrade  his  performance.  The  problems  associated  with  acceleration  forces  acting 
on  pilots  became  even  more  pronounced  with  the  development  of  the  supersonic  high- 
performance  aircraft. 

There  is  a  wide  spectrum  of  acceleration  environments  that  may  be  encountered 
in  a  flying  aircraft;  the  direction  of  the  acting  acceleration  constitute  its  major 
distinctive  characteristic.  A  comparative  nomenclature  for  the  various  conventions 
used  to  describe  the  acceleration  environment,  largely  referred  to  as  G,  may  be 
found  in  the  literature  (Roth  et  al.  [2]).  The  subject  of  the  current  study  was  the 
special  linear  acceleration  which  is  known  as  positive  acceleration,  i.e.,  an  accelera¬ 
tion  acting  in  the  head-to-foot  direction.  Often,  it  is  referred  to  as  +G  or  vertical 

z 

accelerat ion. 

A  vast  amount  of  research  has  dealt  with  the  effects  of  linear  accelerations 
on  pilots.  These  studies  can  be  divided  into  two  categories:  Elementary  research 
pertaining  to  the  psycho-physiological  response  of  the  human  under  G-stress,  and 
research  dealing  directly  with  the  effects  of  sustained  acceleration  on  human  perfor¬ 
mance,  as  expressed  explicitly  by  degraded  tracking  capabilities  and  loss  of  flight 
control.  A  survey  of  these  research  efforts  follows  in  a  subsequent  section. 

4- 

According  to  Brown  [1),  pilots  In  the  Schneider  Trophy  Race  of  1930  were  the  first 
to  have  experienced  the  visual  effects  of  sustained  acceleration. 
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Hies.'  past  rese. ireli  efforts  were  mostlv  empirical.  Some  descriptive  models  of 
human  response  to  f-stress  have  been  postulated,  usually  via  experimental  data  reduc¬ 
tion  uni  statistical  analysis  techniques.  In  order  to  complement  these  research 
eftorts,  and  to  evaluate  pilot  performance  under  G-stress  analytically,  a  normative 
pertormance  model  has  become  a  developmental  goal.  That  is,  an  analytical  model, 
capable  ot  predicting  human  performance  decrements  under  G-stress  is  sought.  This 
tvpc  ot  a  prescriptive  model  can  be  utilized  in  a  variety  of  applications,  such  as 
the  optimization  of  fighter-cockpit  design,  or  the  definition  of  aircraft  state 
augmentat  ion  requirements. 

Acc ’rdinglv,  a  research  program  supported  bv  the  USA F  and  performed  concurrently 
at  tin-  .  . •  :  .rce  Aerospace  Medical  Kern-arch  Laboratory  and  at  CYBKKLAB  the 

1'niverstty  of  Connecticut,  was  undertaken.  Through  this  combined  experimental /analyt ical 
program,  we  have  sought  to  investigate  and  model  the  following  issues: 

i  '  The  effects  of  G-stress  on  human  per  fov.iuncv  in  a  longitudinal  air-to- 
air  compensatory  fracking  task.  Specifically,  experiments  were  carried 
out  with  closed-loop  sustained  accelerations  where  the  subject's  pitch 
rate  commands  were  also  the  driving  signals  to  the  centrifuge  on  which 
the  positive  acceleration  was  simulated. 

I  >  _fhc _ .-Meets  ot  G-stress  on  pilot  pertormance  of  a  secondary  task.  Along 

with  the  longitudinal  air-to-air  tracking  under  G-stress,  the  subjects 
performed  a  lateral  tracking  task  as  well.  This  task  represented  essen¬ 
tially  a  side-task  as  there  were  no  target  aircraft  motions  in  the 
horizontal  plane;  the  tracking  error  was  induced  solely  bv  the  pilot’s 
involuntary  roll  commands.  The  experimental  results  obtained  at  AMR1.  show 
a  significant  degradation  in  performing  this  secondary  task  when  under 
G-stress.  Our  analytical  efforts  duplicate  these  results,  thus  validating 
the  pertormance  model  that  has  been  developed. 
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(3)  Changes  in  pilot  describing  function  when  operating  under  G-stress.  In 
the  second  year  of  the  project  a  new  set  of  experiments  was  devised. 

These  experiments  involved  steady-state  tracking  under  G-stress  and  they 
facilitated  frequency  domain  analysis  of  human  response.  Both  the  experi¬ 
mental  and  the  analytical  phase  of  this  work  explore  the  changes  in  the 
pilot  transfer  function  (between  the  observed  tracking  error  and  the 
control  stick  input)  under  sustained  G-stress. 

(4)  Changes  in  tracking  scores  and  pilot's  input-uncorrelated  control 

( remnant 1 .  One  facet  of  performance-degradation  under  G-stress  is  higher 
RMS  values  of  tracking  error  and  of  control  remnant.  It  is  shown  by  model- 
data  comparison  that  the  performance-model  parameters  are  perturbed  from 
their  nominal  values  when  G-stress  is  introduced. 

(5)  Model  parameters.  A  comprehens i ve  identification  program  was  developed  at 
CYBERLAB  and  greatly  facilitated  the  development  of  the  G-model.  The 
algorithm  was  implemented  in  software  on  our  PDP  11/60  computer  and 
provided  a  mechanized,  efficient  and  objective  means  for  identifying  the 
values  of  the  model  parameters. 

Prior  to  addressing  these  sub-tasks,  a  literature  survey  pertinent  to  the  effects 
of  sustained  accelerations  on  human  beings  was  undertaken. 

1 . 2  Review  of  Physiological  and  Tracking  Performance  under  Sustained  Accelerati ons . 
1.2.1  Cardiovascular  Response:  Perhaps  the  most  fundamental  of  a  person's 
physiological  response  to  sustained  positive  G^-stress  is  the  increase  in  the 
hydrostatic  pressure  gradient  in  the  vascular  system.  The  blood,  being  the  most 
mobile  tissue  in  the  human  body,  is  shifted  from  the  upper  part:;  of  the  body, 
notably  the  head,  and  is  concentrated  in  the  abdomen  and  the  lower  extremities. 
These  observations  were  recognized  as  a  serious  problem  bv  Code  et  al  I  i  , 
and  by  Wood  et  al.  [4],  and  further  confirmed  by  Lindberg  and  Wood  l^l>  and 
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Fraser  [6].  These  articles  report  a  consistent  pattern  of  changes  to  the  cardio¬ 
vascular  system  when  subjected  to  G-stress.  According  to  these  studies,  there 
are  two  distinct  periods  in  the  card iovascular  response  to  positive  acceleration 
with  onset  rate  of  about  1  G/Sec:  (1)  a  period  of  progressive  failure,  and  (2) 
a  subsequent  period  of  compensation.  During  the  cardiovascular  failure  the 
following  phenomena  occur: 

1.  Decrease  in  arterial  pressure  at  head  level; 

2.  Increase  in  arterial  and  venous  pressures  at  the  lower  extremities; 

3.  Increase  in  intra-abdominal  and  intrarectal  pressures  (approximately 
1  inear  with  G^)  ; 

4.  Increase  in  heart  rate  (this  observation  was  also  ascertained  by  Little 

et  al.  [7],  and  by  Smedal  et  al .  [8],  although  the  latter  considered 

mainly  linear  acceleration  in  the  transverse  direction); 

5.  Decrease  in  ear-blood  volume; 

6.  Decrease  in  the  arterial  pressure  pulse  amplitude  at  ear  level; 

7.  Severe  interference  with  blood-tissue  oxygenation,  vision,  cerebral 
activity  and  nervous  system  functions. 

Following  this  period  (~  7  sec)  is  a  subsequent  period  of  compensations. 
Pressor  reflex  triggered  by  the  fall  in  arterial  pressure  in  the  carotid  sinus 
initiates  partial  recovery;  it  produces  hypertension  at  the  heart  level  that 
brings  the  hydrostatic  pressure  distribution  in  the  vascular  system  close  to 
normal  levels.  These  dramatic  changes  in  the  cardiovascular  system  profoundly 
affect  other  physiological  systems,  as  described  below: 

1.2.2  Vision:  Human  visual  capabilities  are  extremely  succeptible  to  variations 
in  the  eye  blood  supply.  Degradation  of  vision  is  further  agravated  by  the 
additional  reduction  in  blood  circulation  in  the  eyes  due  to  changes  in  the 
intraocular  pressure  gradient.  The  gross  effects  of  positive  sustained  accelera¬ 


tion  on  vision  are  as  follows  [9j: 
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L.  An  immediate  consequence  of  decrement  in  the  eye  blood  supply  is 

retinal  hypoxia.  Photographs  of  the  retina  under  sustained  positive 
acceleration  can  be  found  in  an  article  by  Newson  and  Leverett  [10); 

2.  Cross  limitations  in  visual  fields:  tunnel  vision  (loss  of  peripheral 
vision),  greyout,  blackout.  This,  too,  is  a  direct  result  of  the 
reduction  of  blood  circulation  in  the  eye.  This  phenomenon  has  been 
observed  and  reported  by  many  researchers  (White  [11],  York  [12] 
et  al  . ) ; 

1.  Impairment  of  visual  accuitv,  as  reported  by  Frankenhaeuser  [13],  White 
and  lorve  [14],  and  others; 

4.  Increase  in  brightness  discrimination  thresholds,  reported  by  Braunstein 
and  White  [lb]; 

5.  It  has  been  observed  by  Smedal  et  al.  [8]  that  pupillary  dilation  occurs 
with  the  loss  of  peripheral  vision. 

6.  The  oculomotor  mechanism  is  affected.  According  to  White  and  Monty 
[16]  limited  ocular  mobility,  under  sustained  positive  acceleration, 
has  been  observed  due  to  mechanical  changes  in  eye  orbital  tissues. 

These  limitations  can  be  overcome  (with  considerable  effort),  but  may 
result  in  ataxic  eye  movements. 

7.  Smedal  et  al.  [8]  furnished  photographs  of  the  cornea,  under  transverse 
acceleration.  It  is  evident  that  it  is  mechanically  deformed,  resulting 
perhaps  in  blurred  vision. 

1.2.3  Pulmonary-Respiratory  Response:  Physiological  phenomena  related  to  pro¬ 
longed  accelerations  occur  in  the  pulmonary  and  respiratory  systems.  Along  with 
the  increase  in  the  pressure  differences  in  the  vascular  system,  there  is  an 
independent  increase  in  the  arterial  and  venous  pressure  gradients  along  the 
lungs'  base-to-apex  axis.  These  findings  are  mentioned  by  Lindberg  and  Wood  [5] 
and  Fraser  [6].  This  increase  is  followed  by  hypoxemia,  which  increases  in 
severity  with  repeated  exposures  to  vertical  G,  as  reported  by  Barr  [17]. 

Lindberg  and  Wood  [5]  observed  that  displaced  internal  organs,  under  sustained 
positive  G,  pull  the  diaphragm  downward,  causing  difficulties  in  breathing.  York 
et  al.  [12]  and  Smedal  et  al.  [8]  reported  evident  cases  of  dyspnea  (shortness 
ol  breath)  under  transverse  accelerations  (+G  ) .  Smedal  et  al .  [8]  also 
reported  reduction  in  the  vital  capacity,  increase  in  respiratory  rate,  diminished 
alveolar  ventilation  (hypoxia)  and  hypercapnia. 

1.2.4  Psychomotor  Performance:  It  is  quite  obvious  that  under  sustained 
acceleration,  motor  performance  (as  expressed  in  manual  movements  and  manipula- 
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tion  tasks)  would  be  disturbed.  Indeed,  Brown  and  I.echner,  in  their  survey  oi 


research,  report  changes  in  the  forces  that  both  arms  can  exert  when  under  G- 
stress  (1).  Canfield  et  al.  [18]  conducted  experiments  proving  that  movement 
time  and  reaching  errors  increase  with  increase  in  G.  Similar  results  were 
shown  by  Cohen  [19,20],  Kaehler  and  Meehan  [21],  Little  et  al .  [7],  and 

Frankenhaeuser  [13]. 

An  increase  in  human  reaction  time  to  various  stimuli  comprise  another 
symptom  of  sustained  acceleration.  These  findings  are  reported  by  several  re¬ 
searchers  (Frankenhaeuser  [13],  Canfield  et  al.  [22,23],  Kaehler  and  Meehan 
[21].) 

1.2.5  Other:  The  following  are  among  some  other  side  effects  of  positive  G: 

1.  Degradation  in  intellectual  performance,  resulting  probably  from 
acceleration  effects  on  the  central  neural  system  (Grether  [9], 
Frankenhaeuser  [13],  Ross  and  Chambers  [24],  Miller  et  al .  [25]). 

2.  Variations  in  time  perception  with  different  levels  of  G  (Frankenhaeuser 
[13]). 

3.  Fffects  on  the  vestibular,  kinesthetic  and  proprioceptive  senses 
(Brown  and  Lechner  [1]).  These  are  attributed  mainly  to  the  centrifu¬ 
gation  side  effects,  and  are  usually  not  present  during  linear  accelera¬ 
tion  in  real  flight  situations. 

1.2.6  Tracking  Performance 

It  is  our  conviction  that  some,  if  not  most,  of  the  psycho-phys iological 
responses  to  sustained  positive  acceleration  lead  ultimately  to  variations  in 
tracking  task  performance.  It  was  necessary,  therefore,  to  review  the  findings 
related  to  this  problem. 

First  reports  concerning  tracking  performance  and  flight  control  did  not 
appear,  as  noted  by  Grether  [9],  until  1958.  Since  then,  research  in  this  area 
has  been  extensive,  although  no  modeling  efforts  using  modern  control  theory 
have  ever  been  attempted.  The  majority  of  these  studies  has  been  oriented  toward 
transverse  (-l-G^)  direction.  Some  studies,  however,  do  deal  with  tracking  perfor¬ 
mance  (not  necessarily  air-to-air  tracking)  under  +G  stress. 
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Brown  and  Collins  [26]  and  Brown  [27]  pioneered  in  closed-loop  centri¬ 
fuge  simulation  of  air-to-air  tracking  tasks.  They  found  a  significant  increase 
in  what  they  called  "average  tracking  scores,"  mostly  in  the  longitudinal  axis, 
when  performing  under  vertical  acceleration  stress.  Piranian  [28]  and  Smiles 
[29]  also  reported  that  sustained  +G  up  to  6G  had  appreciable  effects,  degrading 
air-to-air  tracking  performance.  Similar  results,  for  an  open-loop  centrifuga¬ 
tion  with  various  levels  of  +0^,  in  pursuit  tracking  experiments  were  reported 
by  Fletcher  et  al.  [30].  Kaehler  [31],  Little  et  al.  [7],  and  Middleton  et  al. 
[32]  conducted  various  tracking  experiments  and,  although  they  considered  trans¬ 
verse  accelerations  only,  their  results  generally  agree  with  the  +G experiment 
results . 

Creer  et  al.  [33]  conducted  tracking  experiments  that  simulated  a  re-entry 
vehicle  piloting  task.  One  of  their  findings  was  that  an  increase  in  tracking 
deficiency  occurred  when  performing  under  +G ^  acceleration,  vs.  +G^  or  -G^ 
stress  conditions.  In  addition,  they  observed  a  marked  deterioration  in  pilot 
tracking  performance  during  the  onset  of  the  acceleration,  i.e.,  during  active 
G.  They  explained  this  by  "vertigo  sensations  caused  by  the  angular  rotations 
of  the  centrifuge  gondola  as  the  centrifuge  was  brought  up  to  the  desired 
operating  speed." 

There  is  a  general  agreement  among  the  various  reports,  relating  to  the 
effects  of  acceleration  feedback  on  tracking  performance.  This  type  of 
informational  feedback  is  usually  categorized  as  motion  cues.  Chambers  [34] 
reported  a  noted  improvement  in  performance  when  the  centrifuge  was  engaged  in 
a  closed-loop  operation  vs.  static  conditions.  However,  this  was  found  for  low 
levels  of  acceleration  only.  Brown  and  Collins  [26]  reported  pronounced  improve¬ 
ment  in  dynamic  vs.  static  tracking  scores  when  a  coordinated  flight  simulation 
was  employed.  Guercio  and  Wall  [35]  also  emphasize  the  importance  of  congruent 
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motion  in  compensatory  tracking.  They  showed  that  the  presence  of  motion  cues 
resulted  in  lower  tracking  errors. 

1 . 3  Experimental  Phase 

Experiments  investigating  the  effects  of  positive  acceleration  on  pilots  were 
conducted  as  early  as  the  mid  1930's.  These  studies  were  carried  out  in  real  aircraft. 
Continuing  research  interests  motivated  the  construction  of  the  first  human  centrifuge 
systems,  for  aviation  studies,  in  the  U.S.  and  in  Germany.  The  basic  human  centrifuge 
is  a  gimbal-mounted  cab  at  the  end  of  a  rotating  arm.  The  cab  is  free  to  rotate, 
independently  of  the  rotating  arm,  in  two  axes.  It  is  possible,  therefore,  for  a  pilot 
seated  in  the  cab  to  perform  a  task  while  subjected  to  a  wide  variety  of  acceleration 
patterns.  Later  models  of  the  human  centrifuge  were  able  to  replicate  air-to-air 
combat  situations  successfully.  A  detailed  description  of  a  human  centrifuge  can  be 
found  in  the  literature  [4,26] . 

One  way  of  inducing  a  sustained  positive  (z-axis)  acceleration  on  a  centrifugated 
subject  is  to  spin  the  apparatus  around  the  central  axis  of  rotation,  while  tilting 
the  cab  inward.  This  is  shown  in  Figure  1. 


Fig  1  POSITIVE  ACCELERATION  INDUCED  ON  A  CENTRIFUGED  SUBJECT 


13 


In  order  to  achieve  the  desired  (or  commanded)  vertical  acceleration, 
and  the  radial  velocity  a>  are  adjusted  according  to  equations  1  . 1- 1 . 2 

radius  of  rotation: 

2 

a  R  cos  t  +  g  sini  -  (1 

z 


the  angle  » 

K  being  the 

(1.1) 


u~R  sini  +  g  cosu  =0  (1.2) 

This  is  the  technique  employed  in  the  high-G  experiments  that  were  conducted  on  the 
Dynamic  Environment  Simulator  (DES)  facility  at  the  Air  Force  Aerospace  Medical 
search  Laboratory . 

The  AMKL  experiments  simulated  an  air-to-air  combat  scenario.  The  centrifuged 
subjects  were  instructed  to  track  a  stvlized  target  aircraft  image  using  a  cross-hair 
reticle  sight.  The  target  motion  followed  either  a  predetermined  or  a  pseudo-random 
pattern  in  either  the  longitudinal  or  the  lateral  plane.  The  closed-loop  configura¬ 
tion  of  the  tracking  task  is  schematically  depicted  in  Figure  2. 


TARGET  INPUT 


Fiq.  2  CLOSED-LOOP  TRACKING  CONFIGURATION 


In  order  to  segregate  the  effects  that  (I  forces  have  on  the  pilot's  tracking 
performance,  two  experimental  conditions  are  employed.  In  the  first  case,  denoted 
here  as  static-*.,  the  subjects  tracked  the  target  with  the  centrifuge  at  rest.  In 
the  second  experimental  condition,  denoted  here  as  dynamic-!,,  the  f)ES 

was  in  motion.  in  addition,  a  third  condition  was  included  in  one  set  of  experiments 

(Chapter  III).  In  these  experiments  the  cumulative  C-effects  were  investigated.  Thus, 

every  G-stress  run  was  immediately  followed  by  a  static  tracking  session.  These 
were  denoted  as  post-static  runs. 

The  recorded  data  included  the  time  histories  of  tracking  errors,  the  subject- 
commanded  and  attained-G,  his  control  inputs  and  other  .auxiliary  variables.  In  addition, 
physiological  variables  were  recorded,  heart  rate,  systolic  and  diastolic  pressures, 
etc.  fhe  raw  data  analysis  included  ae ross-sub j ect  ensemble  averaging  and  EFT  analysis 
on  stationary  data. 

A  detailed  description  of  the  experimental  program  as  well  as  of  the  data 
analysis  techniques  involved  follows  in  Sections  2.1  and  3.1. 

1 .  A  Mode  l_in^  Cons  i de  rat  Ions 

The  modeling  efforts  in  the  present  study  have  been  based  upon  the  well-known 
Optimal  Control  Model  (OCM)  (K.leinman  et  al.  [36],  Baron  et  al.  [37],  and  Kleinman 
et  al.  [18]).  The  main  contribution  of  the  modeling  phase  was  the  adjustment  and  the 
modification  of  the  OCM  to  Include  the  effects  of  G^-stress  on  pilot  performance.  A 
major  reason  for  using  the  OCM  as  opposed  to  some  other  human  performance  model,  is 
that  the  parameters  of  this  model  are  descriptive  of  human  limitations.  Thus,  they 
have  shown  a  remarkable  invariance  over  different  task  specifications,  input  distur¬ 
bances  or  commands,  and  controlled-vehicle  dynamics.  It  therefore  follows  that 
changes  in  the  "nominal"  OCM  in  a  sustained  acceleration  environment  are  induced 
almost  entirely  bv  the  G  stress. 
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.  .i  Modeling  Approach 

1  shews  the  si  rui'i  u:r  <>t  t  he  OCM.  Tl*t-  veil  i  «  I  e  dynamics  are 
1  -  in  sl.d  •  spare  farm 

lt  >  -  A  >.  ft)  f  b  u(t  >  +  I.  w  (t)  +  K  /.(l  i 

O  C.— O  -o  l '  O  <* 

v'l.)  =  C  *  (t)  +  d  uli)  (  1 .  i) 

where  ••  •  t )  is  usually  a  white-Gaussi.in  disturb, im  e  with  covariance  W  and 

.i  o 

:  a  i  ■_  rmn  i  t  K  vehicle  r  ru  j  t  ct  •  >r  y .  The  displayed  information  y(t) 

onsi-  h  ,  t  t  •  r.;,:  error  eft)  aid  error-rate  eft),  as  wel  1  as  other  auxiliary 

variables  and  their  rate  of  change-. 

The  visual  assumption  in  the  OCM  is  that  the  human  perceives  a  delayed  and 
noisy  replica  of  y(t),  viz., 

_y  (t)  =  _y(t— t)  Vy f t  —  t )  (1-4) 

where  the  white-Gaussian  observation  noise  v^(t)  °f  the  i-th  indicator  has 


:ovarianee 


(t)  =  -Xi _ _  Et2(t)l 

rr  }  f.(t)N(a.) 


In  these  equations  t  is  the  pilot's  lumped  time-delav,  p  is  the  observation 
noise/signal  ratio,  f^(t)  is  the  fraction  of  attention  the  human  is  allocating 
to  each  observation^,  and  N(a^)  is  the  equivalent  "gain"  of  the  visual/indif¬ 
ference  threshold  a^.  The  gain  N(a^)  is  obtained  via  statistical  linearization 
derived  from  Random  Input  Describing  Function  Theory,  and  is  also  a  function  of 
the  mean  and  the  standard  deviation  of  y^(t). 

Assuming  that  the  system  state  estimate,  x^,  is  given,  then  the  human 
develops  an  optimal  control  strategy  by  minimizing  the  quadratic  cost  functional, 


Total  attention  f  is  assumed  =  1.  This  is  subdivided  across  axes  of  control,  and 
further  subdivided  among  instruments  per  axis. 
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OPTIMAL  CONTROL  MODEL  OF  HUMAN  RESPONSE 

DISTURBANCES 


OPTIMAL  CONTROL  MODEL  OF  HUMAN  RESPONSE 


'V 


J(u)  =  £im  ^  E  {  I  (y'  (t)Qir  y(t)  +  gii2(t))  dt 


(1-6) 


where  Q =  diag  (q^)  are  the  relative  weightings  on  the  observations  and  g 
is  the  control-rate  weighting^.  Usually,  the  tracking  error  weighting  co¬ 


efficient,  q^,  is  non-zero,  since  the  objective  is  to  minimize  the  observed 
error.  By  augmenting  the  state  estimate  with  the  control  u,  viz., 


(1.7) 


this  strategy  results  in  the  optimal  control  gains,  L^,  arising  from  the 
pertinent  steady-state  Riccati  equation,  such  that 


u(t)  =  -4  °1  =  -L’*(t) 

if  we  define  l!  =  It.,1  l'  .  T.,1  and  include  motor  noise  in  the  pilot's 
-1  [  N  —opt  N  J 

control,  the  optimal  control,  u(t),  obeys  the  equation 

V  +  u  =  "-opt  4  +  Vu(t)  “  Uc(t)  +  Vu(t)  (1-8) 

or  with  l'  =  ("  L  J  ol 

-  L“°Pt!  J 

ru+u=L'x+v(t)=u(t)+v(t)  (1.9) 

N  -  u  c  u 

The  parameter  iN  can  be  interpreted  as  a  "neuromotor"  time  constant.  Usually, 

t  ,  is  specified  and  g  is  adjusted  accordingly,  as  there  is  a  one-to-one 
N 

correspondence  between  the  two. 

The  motor  noise  v  (t)  is  assumed  to  be  white  and  Gaussian  with  covariance 


We  assume  that  the  subject  seeks  a  control  strategy  that  would  minimize  control 
rate  rather  than  controL. 
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that  scales  with  the  covariance  of  u(t), 


V  (t)  =  Ttp  cov  [u(t)  1 . 
u  u 


(1.10) 


The  coefficient  represents  the  motor-noise/signal  ratio. 

Next,  it  is  necessary  to  examine  the  other  blocks  of  the  human  information 
processing  unit  (figure  2).  That  is,  how  does  the  subject  estimate  the  state 
jt(t-x),  and  how  does  he  predict  £(t)  from  x(t-x).  The  augumented  state 

x'(t)  =  (x'(t),  u(t) ]  obeys  now  the  differential  equation 

—  — o 


x(t)  =  A  x(t)  +  bu^(t)  +  Ew(  t )  +  Fz(t) 


(1.11) 


where 


'  k 
!  ° 

-  !"  “-I  ’  - 

0  :  J 


o  iv1 


"f 

“O 

0  J 


(1.12) 


and  w(t) 


wo(t) 


vu(t) 


with  the  covarience  W  = 


(1.13) 


It  is  not  difficult  to  show,  that 


x(t)  =  e^T.x(t-T)  +  J  e  ^  buc(o)  do 


(1.14) 


InEquation  (1.14) the  subject  will  use  the  estimate  x(t-x)  =  p(t).  This  state 
estimate  is  obtained  from  the  Kalman  filter,  via 


p  (t)  =  Ap(t)  +  bu  (t-x)  +  C(t)  [v  (t)  -  Cp(t)  ] 
x-  —  —  c  ~T) 


(1.15) 
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(1.16) 


where 


p(t)  =  E  |  x(t-x) llp(o) .  o<t| 


The  filter  gain  G(t)  is 


C(t)  =  E(t)C*V  X(t)  ,  Vy(t)  =  diag  [V  (t) ] 


(1.17) 


where  E(t)  is  generated  from  the  solution  of  the  Riccati  equation 


E  =  AE  +  EA'  -  EC ' V  1CE  +  EWE'  +  W^t) 
y  d 


(1.18) 


Three  observations  should  be  made  with  regard  to  Equation  (1.18). 


1.  The  matrix  W^(t)  is  a  "pseudo-noise"  covariance  matrix  and  arises  from 
the  deterministic  signal  z(t).  It  is  chosen  in  such  a  way  that  the 
filter  adapts  to  the  changing  characteristics  of  the  trajectory  z(t). 
At  present, 

W,(t)  -  x  F  z2(t-T)F'  (1.19) 

d  cor  —  — 

where  icor  is  linked  to  the  human  short-term  memory,  or  alternately, 
is  dependent  on  the  z(t)  profile's  bandwidth. 

2.  In  the  absence  of  the  deterministic  target  motion,  W^(t)=0  and  E=0. 
That  is,  the  steady-state  Riccati  equation  is  solved. 


3.  In  the  case  when  the  noisy  disturbance  w^ftJ'O,  the  term  EWE'  becomes 


0 

0 


i  0 


(1.20) 


1.4.2  Model  Application 

Basically,  modeling  efforts  utilizing  the  OCM  can  be  approached  in  two 


ways : 

1.  Time-varying  (nonstationary)  mode 

2.  Steady-state  mode. 


20 


1.4.2  .JL_  Nonstnt  i  onai  y  Mod.  1  • 

In  the  time-varying  case,  the  target’s  profile  is  deterministic, 

i.e.,  z(t)^0  and  w  (t)=0.  This  model  tug  to*  i-r  ique  usin'.  i  invc- ;  vt*s  a 
o 

compel  i  son  of  exper  j  merit  a  1  t '  nv  his-  ■■  u> .  i  1 1  ■>  -  ■■■  .  ■  ,t  •. 

statistics)  with  the  mode  1- ..  rod  it  ted  at.  *>  1  :  /  ;  •  •,  ....ert-i  o 

to  develop  t.  u  process  me-.n  ami  c.ovat  tar.-  ptopag.t  iou  <n;u«t  io;»s  whit, 
at  use  from  the  non-rand--, >•  i  ,  mpor-enc  z.h  i  .  -t  ..  ■  ...  .  * . •  n .  ihese 

equations  have  been  thoroughly  documented  :  e.  g.  I  18  i  ,  [  "1°  .  I  4o  j ; .  and 
the  basic  continuous  time  equations  are  repeated  here  wlturti  proofs  for 
completeness . 

If  one  defines 


x  =  E(x) 
x  -  E(x) 

jg^(t)  =  x^t-r)  -  £(t)  =  filtering  error  ;  =  E(e^)  (1.21) 

^2(t)  =  x(t)  -  x(t)  *  prediction  error  ;  e^  =  Efe^) 

1  =  E(f.) 


then  the  process  mean  and  covariance  are  computed  from  the  following 
time  equations: 


|(t)  =  A  J(t)  +  G  C  e^t)  ;  A  =  A-bl.'  ,G  =  eATC, 


e1(t)  =  A  e^t)  +  F  z(t) 


A  =  A-GC 


l: 


e2(t)  =  /  eA(t  a)  F  z(a)da  ; 


x(t)  -  |(t)  +  eAT  e^t)  +  e2(t) 


i(t)  =  Cx(t) 


(1.22) 

(1.23) 

(1.24) 

(1.25) 

(1.26) 
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The  boundary  conditions  are  appropriately  chosen  at  time  t=T . 


Next,  if  we  let  X,  E  ,  E,,.  X  and  Y  be  the  covariance  matrices  of 


x,  e  ,  e x  and  y  respectively,  and 


M(t)  =  E  (x(t)  -  x(t))  (ej(t)  -  e 


one  can  show  that 


j  (O)'  - 


cross  covariance,  (1.27) 


X(t)  =  X(t)  +  eATM'(t)  +  M  (t)eA  T  +  eAlEL(t)eA  X  +  E2(t)  (1.28) 


E2(t) 


t 

=  f  eA(t-a)  EW(o)E*  +  Wd(a)j  eA  (t‘ 

Jt- T 


(1.29) 


Fn(t)  =  AE^t)  +  E1(t)A  +  GV  (t-T )  G  ’  +  EW(t~x)E’  +  Wd(t-T) 
M(t)  =  AM(t)  +  M(t)A  +  GC  [Ei(t)  -  F.(t)j 


X  =  AX  +  XA'+  Gen' (t)  +  M (t)c'c'  +  GV  (t-T)c' 


(1.30) 


(1.31) 


(1.32) 


Y  =  CXC 


(1.33) 


Notice  that  since  we  assume  in  this  case  w  (t)-0,  the  matrix  W(t)  is  a 

o 

zero  matrix  except  for  the  lower  right  lxl  corner  entry  which  is  vu(r) 

(see  Section  1.4.1). 

I. 4. 2. 2  Steady-State  Mode 

In  the  steady-state  mode  z(t)=0,  and  the  target  input  is  driven 
by  a  stationary  colored  noise.  The  usual  practice  is  to  describe  the 

disturbance  w  (t)  as  a  white-Gaussian  noise,  and  to  augment  the  system 

o 

with  the  dynamics  of  a  noise  shaping  filter  that  charai terizes  the 
driving  noise.  The  stationary  approach  is  most  suitable  for  frequency 
domain  analysis  of  human  response.  The  model  equations  are  represented 
in  the  frequency  domain,  so  that  various  performance  measures  can  be 
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predicted.  The  model  outputs  include  any  of  several  possible  describing 
functions  associated  with  the  subject,  the  vehicle,  or  the  overall 
closed-loop  system.  Also  available  are  power  density  spectra  of  any 
pertinent  system  variable  that  show  both  correlated  (i.e.,  input-related) 
and  remnant  (i.e.,  human- related )  components.  In  addition,  the  perfor¬ 
mance  scores  (RMS  values)  are  predicted.  Among  the  quantities  that  are 
considered  in  the  present  study  are 

1.  Pilot  describing  function, 

2.  Pilot  control  remnant,  and 

3.  Tracking  scores:  RMS  values  of  the  tracking  error  and  control. 
This  modeling  approach  is  most  suitable  for  model  identification 

from  frequency  domain  data.  The  model  is  applied  in  conjunction  with 
f ourier-transformed  experimental  data,  and  measures  such  as  described 
above  can  be  compared. 

The  frequency  domain  equations  are  summarized  here.  A  more  detailed 

theoretical  background  can  be  found  in  the  literature  (e.g.  [41]). 

We  assume  the  model  of  Equation  (1.11)  with  z(t)=0.  All  initial 

conditions  are  ignored  since  steady-state  situation  is  assumed.  The 

major  transfer  function  of  interest  initially  is  that  between  y(s)  and 

u  (s)  : 
c 

U  (s)  =  H1  (s)  jr(s)  (1. 34) 

c  — 

This  function  is  obtained  by  taking  the  Laplace  transform  of  the  equations 
of  Section  1.4.1.  The  following  relationships  result: 

H'(s)  =  [l  +  r(s)]_1L’e(A"Sl)T  (sI-A)_1G  (1.35) 

r(s)  =  l’|  e(A_Sl)T  [(sI-A)'1  -  (Sl-A)_1j  +  (s!-A)_1|b  (1.36) 
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Tlu'  overall  input-output  human  transfer  function  is  thus 

u(s)  =  -(s  +  tn  S  Hi  (s)  y(s)  =  I  '  (s)v(s)  (1.37) 

where  we  note  that  (s4r^  )  is  obtained  I  rom  the  lower  right  lxl 

corner  of  (sI-A)  Next,  we  also  require  the  decomposition  of  all 

closed-loop  signals  jc,  y  and  u  into  their  two  major  components.  The 
first  is  that  part  linearly  correlated  with  the  independent  input 
disturbance  w^(t).  The  second  is  that  part  due  to  the  human-injected 
remain t  y^  and  v^.  Thus,  with  ^  =  [E^,  0),  we  obtain 

u,(s)  -  -li'(s)Cx(s)-  h'(s)v^(s)+  (1.38) 

x(s)  -(si-A+bli  C)  'k  w(sl  +(sf-A+bN,C)ib  w'v  (s)+t  (s)  (1.39) 

—  -o  v  —  -  y  NuJ 

v(s)  -  Cx(s)  (1.90) 


following  these  equations  the  human  transfer  function  of  (1.37)  can  be 
obtained.  This  describing  function  is  of  interest  to  us  since  it  can  be 
measured  experimentally.  It  is  defined,  relative  to  the  disturbance 

w  (t),  as  a  ratio  of  the  transfer  function  between  u  and  noise  w  to 

o  o 

the  transfer  function  between  the  observed  tracking  error  e  and  noise 

w  .  Thus 
o 


u(s_) 

e(s) 


u(s)/w  (s) 
o 

e  ( s )  /w  ( s ) 
o 


(1.91) 


To  compute  this  human  describing  function  we  use 


w"  7s)  -  j^I-A+bHV^ 

(1.92) 

:7,  - 

o  Ju 


It  is  recognized  that^Vy(t)  does  not  strictly  have  a  Laplace  transform.  This, 
however,  is  not  a  problem  since  we  are  interested  in  the  uncorrelated  spectra. 
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The  subscripts  e  and  u  in  (1.42)  indicate  that  the  tracking  crror- 
and  control-related  components  of  the  describing  function  are  considered. 

Finally  we  require  power  spectral  densities  of  the  various  variables 
in  the  loop.  These  are  obtained  from  liquations  (1.38-!. 40).  lor  anv 
signal  the  spectrum  is  decomposed  into  input  and  remnant  related  com¬ 
ponents,  (superscripts  c  and  ue  respectively)  and  the  spectral  densities 
ij)  (•)  are  obtained,  viz., 


f  ■  ‘  J  !*;,<miv„ 


(1.43) 


=  i!Bk(-i->lWo 


(1.44) 


where, 


k-th  state  component 

ki-th  element  of  (sI-A+bH'c)  ^bH 

,  i  -1 

k-th  component  of  (sI-A+bH  C)  Jjt 


N  =  number  of  observed  variables  (dimension  of  y(t)) 

y 

i  _i 

6^  =  k-th  element  of  (si-A+bH  C)  E 

c 

Similar  expression  are  easily  obtained  for  |C  (to)  ^uc(w)  (j)  (w)  and 
uc  uc  ^k'  ^k  u 

|  (w)  The  spectral  density  |  (w)  is  of  special  interest  in  this 

u  u 

study  as  it  represents  the  human  input-uncorrelated  control.  It  will 
be  shown  in  chapter  Ilf  that  the  magnitude  of  jji^C  ('“)  in  G-stress 
conditions  is  significantly  larger  then  it  is  in  no-G  situations.  This 
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is  an  excellent  indicator  of  the  human's  motor  randomness  which  increases 


in  the  presence  of  sustained  linear  accelerations. 

1.4.3  Model  Development 

The  baseline  model  adopted  in  this  work  was,  as  noted  previously,  the 
Optimal  Control  Model  of  human  performance.  The  approach  taken  in  the  G-§tress 
modeling  process  was  to  develop  a  "G  submodel"  that  can  be  embedded  in  the 
general  structure  of  the  OCM.  The  steps  taken  in  this  direction  included  a 
development  of  technology  to  isolate  and  identify  the  effects  on  tracking 
error  and  pilot  control  of  target  motion  versus  those  associated  with  G-stress. 
The  effects  of  target  motion  were  studies  via  static  (no  G)  experiments.  If 
a  pilot  model  is  first  tuned  accurately  to  this  case,  then  our  hypothesis  was 
that  any  data-vs-model  mismatch  in  the  dynamic-G  tracking  case  is  attributable 
to  the  effects  of  the  G-stress.  The  subsequent  identification  process  to 
model  the  mismatch  consisted  of  four  steps: 

1.  Sensitivity  analysis:  A  comprehensive  sensitivity  analysis  of 
model  to  OCM  parameter  variations  was  undertaken. 

2.  Development  of  generic  structural  formulae  that  relate  pilot  para¬ 
meters  to  instantaneous  and/or  antecedent  load  factor  (G  and  JG) ; 

3.  a.  Adjustment  (tuning)  of  the  free  parameters  within  the  OCM: 

t,  tn,  a.,  py.  and  P^; 

3.  b.  Tuning  of  the  free  parameters  within  the  structural  submodels  so 

as  to  best  match  given  data; 

4.  An  identification  program  that  performed  the  final  fine-tuning 
adjustments  to  the  OCM  and  the  G  submodel  parameters. 

In  addition  to  obtaining  a  performance  model,  this  modeling  approach  provides 
a  potentially  powerful  method  of  obtaining  physiological  correlates  of  perfor¬ 
mance  under  acceleration  stress.  We  believe,  therefore,  that  an  ultimate 
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modeling  hierarch  which  defines  in  succession  the  functional  relationships 
between 

1.  Stress  and  physiology 

2.  Physiology  and  control  parameters 

3.  Control  parameters  and  tracking  performance  (via  OCM) 

will  provide  maximal  flexibility  for  future  work.  This  hierarchy  is  shown 
schematically  in  Figure  4.  It  should  be  noted,  however,  that  even  a  failure 
to  obtain  the  submodels  of  levels  1  and/or  2  does  not  hinder  the  general  success 
of  such  approach.  This  "robustness"  of  the  technique  is  significant  since 
research  in  the  area  of  modeling  physiological  factors  has  traditionally  been 
painfully  slow  and,  to  date,  not  very  fruitful.  It  is  always  possible  to 
consolidate  the  first  two  levels  into  a  single  augmented  G-stress-control 
parameters  model,  circumventing  the  physiological  variables  altogether. 

In  the  following  chapter,  the  experimental  and  modeling  phases  of  this 
project  are  detailed. 
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HUMAN  OPERATOR  MODELING  HIERARCHY 


II.  REPETITIVE  -  G  —  NONSTATIONARY  TARGET  PROF II  .K 

In  the  first  year  of  the  program,  our  efforts  were  directed  towards  analyzing 
the  Repetitive-G  ensemble  data,  followed  by  the  development  of  an  analytical  model. 

The  Repetitive-G  experiments  were  performed  on  the  DES  facility  at  AMRL,  and  pro¬ 
vided  the  data  base  for  our  nonstationary,  time-domain,  G-Stress  performance  model. 

2.1  Experimental  Program 

Repeated  trials  of  a  compensatory  tracking  task  were  performed  on  the  DES  at 
AKKI..  The  subjects  were  seated  in  the  centrifuge  cab  and  were  instructed  to  track  a 
simulated,  triangular-shaped  target  aircraft  image  that  was  displayed  on  a  graphics 
screen.  The  target  motion  followed  a  predetermined  pattern  in  the  longitudinal  plane. 
There  was  no  target  input  in  the  lateral  plane,  but  the  operator  had  to  correct  hori¬ 
zontal  tracking  errors  that  were  induced  bv  his  own  inherent  motor  randomness.  Figure 
5  shows  the  structure  of  the  visual/motion  loops.  G  (s)  ar>d  G^fs)  represent  the 
controlled  vehicle  dynamics  in  the  pitch  and  the  roll  axes  respectively;  these  are 
discussed  later. 


Fiq.  5  CLOSED-LOOP  DYNAMIC  TRACKING  CONFIGURATION 
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.  L_.  1_  Target  Profiles 

The  typical  target  maneuver  in  the  longitudinal  (pitch)  plane  is  represented 
in  terms  ot  the  induced  vertical  acceleration  G^,(t).  Three  different  target  tra¬ 
jectories,  employing  trapezoidaly-shaped  acceleration  profiles  were  used; 

(a)  Mo  Peaks:  Shown  in  Figure  b;  runtime  =  103.34  seconds 

(b)  One  Peak:  Shown  in  Figure  7;  runtime  =  121.33  seconds 

(e)  Five  Peaks:  Shown  in  Figure  8;  runtime  =  260.  seconds. 

The  target  profiles  are  characterized  by  the  following  features: 

I)  The  peak  levels  are  7C  and  are  sustained  for  a  period  of  10  seconds. 

2’  There  is  a  10  second  rest  period  (1.5G)  between  two  successive  peaks. 

‘3>  The  rate  of  onset  of  G-j.(t)  from  the  1.5  G  level  to  the  7G  peak  level  is 
0.75  G/Sec. 

4)  There  is  a  90  second  final  period  of  a  constant  4G  level  plateau  for  each 
of  the  three  profiles.  The  subjects'  performance  was  scored  during  these 
periods  only. 

The  profile  Gj(t)  can  be  interpreted  as  "commanded-G"  input  to  the  human- 
tracker.  Tf  the  human  operator  (HO)  were  to  track  with  zero  error,  G-p(t)  would 
equal  the  acceleration  forces  he  would  experience.  Acturally,  the  subject's 
"attained  vertical  acceleration,"  G^(t),  would  be  a  filtered  version  of  Gf(t) 
because  of  the  lag  inherent  in  the  HO-System-DRS  dynamics.  Since  the  aircraft's 
attained  G  level  is  directly  proportional  to  its  pitch  rate,  one  may  envision  the 
simulated  scenario  as  a  pursuit-evasion  task  in  a  circular  path  with  a  changing 
radius.  The  general  equation  that  relates  the  aircraft  pitch  rate  to  its  attained 
G  is  given  by 
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(2.1) 


G(t) 


YSL (O 

g 


+  1 


where  V  is  the  aircraft's  velocity  (assumed  constant)  q(t)  is  its  pitch  rate,  and 
g  =  )2.2  ft/sec.  The  1  accounts  for  the  (fixed)  acceleration  of  gravity. 

2.1.2  Experimental  Conditions 

The  major  goal  in  the  design  of  these  experiments  was  to  detect  any  perfor¬ 
mance  degradation  under  sustained  G-stress.  Therefore,  two  experimental  condition 
had  to  be  employed. 

1,  Static-G:  The  subjects  tracked  the  target  with  the  centrifuge  at  rest. 

2.  Dynamic -G :  The  l)KS  was  engaged  in  a  closed-loop  mode.  The  subjects, 

responding  to  a  pitch  error  stimulus,  generated  a  command  control  input  in 

order  to  reduce  the  tracking  error.  By  this  action,  they  induced  a  positive, 

time  varying  G  stress  upon  themselves  by  increasing  or  by  decreasing  the 
z 

angular  velocity  of  the  centrifuge  (see  also  Figure  5). 

Note  that  the  roll  error  existed  only  in  the  visual  loop  and  no  off-normal 
(Gy)  accelerations  were  generated  by  the  pilot's  roll  commands. 

As  a  consequence,  a  3x2  factorial  design  was  utilized.  Each  of  the  subjects 
tracked  three  different  trajectories,  and  each  trajectory  was  tracked  both  under 
conditions  of  "Static-G"  and  "i)ynamic-G." 

2.1.3  Subjects 


The  AMRL  hazardous  duty  panel  provided  the  personnel  for  these  experiments. 
Eight  male  subjects  were  trained,  and  each  provided  approximately  six  runs  per 
target  trajectory  per  experimental  condition. 
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2.1. 4  Aircraft  Longitudinal  (pitch)  Dynamics: 

The  basic  set  of  longitudinal  equations  being  used  for  the.  attacker's  air¬ 
craft  is  the  short-period  dynamics.  They  are  the  pertubation  equations  written 

about  straight  and  level  flight,  0  =0. 

o 

a  =  q.  +  Z  a  +  z,  6 .  (2.2) 

A  A  a  A  6  A  ' 


q  =  M  0t.  +  M*  a.  +  M  q,  +  Mx  6. 
A  a  A  aA  q  mA  6  A 


(2.3) 


6  A  =  qA 


(2.4) 


Where  9  is  the  attacker's  pitch  angle,  q  is  the  pitch  rate,  a  is  the  angle  of 
A  A  A 

attack  and  6^  is  the  elevator  defection  (all  pertubation  values).  The  normalized 

stability  derivatives,  Z  ,  Zx,  M  ,  M*  ,  M  and  Mx  are  generally  functions  of  the 

a  6  a  a  q  6 

nominal  angle  of  attack  and  dynamic  pressure.  The  range  of  validity  of 
Equations  (2.2)— (2.3)  with  constant  parameters  may  be  small,  contributing  to 
modeling  errors  in  an  air-to-air  problem.  Also,  the  choice  of  equilibrium  condi¬ 
tion  may  vary. 

In  Equation  (2.2)  the  parameter  Z^  is  usually  small  and  sometimes  neglected. 

The  same  is  sometimes  true  for  M-  .  The  analysis  that  follows  will  assume 

a 

M-  =  Zx  =  0  (2.5) 

a  o 

as  this  corresponds  to  the  assumptions  made  in  the  simulations.  It  should  be 
noted  that  this  does  not  compromise  the  results;  non-zero  choice  presents  no 
analytic  difficulties. 
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The  transfer  functions  q  /6.  and  a  /&  are,  with  Equation  (2.5) 

A  A  A  A 


q.  Mrs  -  MXZ 
mA  _  6 _ 6  a 

6  A(s) 

A 


A. 

A(s) 


(2.7) 


where 


A(s) 


2 

s 


(Z  +  M  )  s  +  M  Z 
a  q  q  ot 


M 

a 


(2.8) 


Acceleration  Forces 

The  acceleration  normal  to  the  flight  path  due  to  an  incremental  flight- 
path  angle  is 


G  =  V  yA 
Z  A 


where  V  =  aircraft  velocity  (assumed  constant)  and 
angle.  Noting  that 

“a  +  V*  ■  \ 


(2.9) 

(incremental)  flight  path 

(2.10) 


or  equivalently 


and  using  Equations  (2.2)  and  (2.11)  one  obtains 

2 

G  =  -  V  Z  a  ft/sec 
z  a  A 

Thus,  the  transfer  function  is 


(2.11) 


(2.12) 


G 

z 


Z  Mx 

a  6 

A(s) 


(2.13) 
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From  the  transfer  function  parameters  given  in  the  HAC  study,  A(s)  =  s  +9.9s+38 
and 

M,  =  11,  Z  =  -25/11,  V  =  1000  ft/sec 

0  ct 

M  =  -7.63,  M  =  -20.66 
q  a 

Geometric  Considerations 

Figure  9  shows  the  geometry  of  air-to-air  tracking  in  the  longitudinal 
plane.  The  variables  with  subscripts  A  relate  to  the  attacker  and  those  with  T 
to  the  target  aircraft.  There  are  no  gunsight  dynamics,  i.e.  the  sight  is  fixed 
and  aligned  with  the  aircraft  body  axis.  For  the  case  of  gunsight  lead,  etc., 
see  Harvey  and  Dillow  [42]. 


r  =  aspect  angle  of  target  =  0j  -  Ej 
D  =  distance  between  the  two  aircraft  =  1000  ft. 


Fig.  9  LONGITUDINAL  TRACKING  GEOMETRY 
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Assuming  V  ~  V  =  V  and  small  (relative)  angles  it  has  been  shown  [42]  that 

A  i 

V  D  (YT  ■  V  (2-14) 

for  constant  V  and  D.  Since  the  explicit  second  order  response  of  the  target  is 
not  important  for  modeling  the  attacker,  we  conveniently  assume 

qT(t)  =  z(t)  (2.15) 

where  z(t)  is  a  function  of  target  profile  as  shown  in  Figures  6-8  and  described 
in  Section  2.1.1.  By  referring  to  Equation  (2.1),  one  may  also  recognize  z(t)  as 

z(t)  =  *  GT(t)  (2.16) 

2.1.5  Aircraft  Lateral  (roll)  Dynamics 

The  lateral  tracking  task  represents  essentially  a  side-task  as  there  are 
no  target  aircraft  motions  in  the  horizontal  plane.  In  order  to  model  the  lateral 
task,  the  dynamics  between  aileron  deflection,  6^,  and  tracking  error,  e^,  must 
be  defined.  The  implementation  of  the  lateral  tracking  task  assumes  the  following: 

1.  Roll  angle,  <J>,  is  small  so  that  cos<}>~l  and  sin<{)~4).  This  enables  us  to  treat 
the  lateral  and  longitudinal  modes  as  uncoupled.  Moreover,  the  longitudinal 
tracking  task  is  not  dependent  on  4>  in  this  case. 

2.  The  attacker  angle-of-attack  a and  sideslip  angle  8^  are  negligible.  Thus, 
the  attacker  velocity  vector  is  always  aligned  with  the  body  axis.  This 
greatly  simplifies  the  system  representation  in  the  horizontal  plane. 

3.  All  turns  are  coordinated,  i.e.  the  G  vector  is  aligned  with  the  aircraft's 
z-body  axis. 
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if  we  consider  y  =  constant  then  Equation  (2.23)  represents  a  dynamic  lineariza¬ 
tion.  In  the  present  problem  v  is  a  random  variable,  with  a  mean  and  variance 

A 

that  are  functions  of  time.  It  is  more  convenient  to  use  y  in  Equation  (2.23), 
since  this  will  be  a  "cleaner"  signal.  The  approximation  is  valid  to  first-order 
as  (Yt  -  Y  )  0  =  0.  Thus, 


=  (yt  +  y) 


(2.24) 


The  next  equation  specifies  the  relationship  between  heading  error  v  and 
lateral  tracking  error  e  (t).  It  is  easy  to  show  that 


e 

r 


(2.23) 


where  D  =  distance  between  the  two  aircraft. 

The  final  equation  reflects  the  aircraft  roll  axis  dynamics  which  is  essen¬ 
tially  a  roll-rate  command  system 


XA  $  +  $  =  6  (2.26) 

<p  a 

In  the  present  simulations  —  =  5. 55  ,  D  =  1000  ft  and  V  =  1000  ft/sec. 

T, 

<P 

Equations  (2 . 24)- (2 . 26)  show  that  the  response  between  6a  and  er(t)  is  like 
3 

1/s  .  This  is  a  very  difficult  system  for  a  human  to  control,  depending  on  the 

visual  information.  The  block  diagram  of  this  system  is  shown  in  Figures  10. 

.  2 
Note  that  Vy^  +  g  is  the  target  normal  acceleration  in  ft/sec  .  For  convenience, 

this  open-loop  gain  has  been  made  into  an  equivalent  gain  K(t)  at  the  input. 

With  these  equations,  the  system  description  in  both  the  longitudinal  (G^(s))  and 

the  lateral  (Gr(s))  axes  is  complete. 


38 


Fiq.  10  A  4th  ORDER  LATERAL  TRACKING  MODEL 


2.1.6  Data  Acquisition 

The  data  supplied  by  AMRL  consisted  of  the  time  histories  of 

1.  Longitudinal  (pitch)  axis  tracking  error 

2.  Lateral  (roll)  axis  tracking  error 

3.  Commanded  G  (input  to  DES) 

z 

4.  Attained  G 

z 

5.  Pilot  longitudinal-axis  stick  input 

6.  Pilot  lateral-axis  stick  input. 

These  were  provided  for  the  0-peak,  1-peak  and  5-peak  trajectories.  Two 
magnetic  tapes  were  provided,  one  for  each  experimental  condition  (Static-G  and 
Dynamic-G) . 

The  sampling  interval  in  the  real-time  operation  of  the  formal  experimental 
runs  was  0.06  seconds.  The  data  were  collected  every  other  time  step,  providing 
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a  sample  of  each  of  the  six  variables  every  time  frame  of  0.12  seconds.  This 
produced,  for  every  recorded  variable,  the  following: 

1.  862  datum  points  in  the  0-peak  profile. 

2.  1012  datum  points  in  the  1-peak  profile. 

).  2168  datum  points  in  the  5-peak  profile. 

Tiie  total  number  of  complete  runs  that  were  suitable  for  reduction  and  analysis 
was : 

-  Static  Condition:  0-peak  ,  N=41. 

1-peak  ,  N=38. 

5-peak  ,  N=32. 

-  Dynamic  Condition:  0-peak  ,  N=37. 

1-peak  ,  N=36. 

5-peak  ,  N=24. 

2.2  Ensemble  Data  Analysis 

The  raw  data  obtained  on  the  DES,  in  conjunction  with  the  repetitive-G  tracking 
experiments,  have  been  analyzed  via  aggregation  using  an  ensemble  averaging  method. 
This  has  been  done  for  subsequent  model-data  detailed  comparisons,  and  for  covariance 
modeling  modifications. 

An  extensive  set  of  software  was  developed  for  processing  the  single-run  data. 
Tape  conversion  programs  were  written  on  UConn  Computer  Center's  IBM  360/65  to  read 
the  AMRL  tapes  and  to  rewrite  them  in  a  form  compatible  with  our  laboratory's  DEC 
PDP-11  system.  The  PDP-11  system  was  used  for  the  averaging  process,  as  it  enabled 
a  hands-on,  interactive  analysis.  Accordingly,  programs  were  written  for  averaging 
the  data  in  a  serial,  or  sequential,  manner.  In  this  way  the  editing,  addition,  or 
removal  of  a  single  run  could  be  done  iteratively,  without  having  to  re-average  the 
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entire  ensemble.  This  is  an  extremely  useful  approach  for  treating  outlyer  runs 


not  satisfying  a  t-test  criterion.  Specifically,  the  following  equations  were 
employed  to  reduce  the  empirical  data  process: 

1.  The  ensemble  mean  x^+j  (i)  of  the  i-th  point  of  k+1  runs  is  computed  from 

\(i)  by 

Xk+l(1)  =  (\(i)‘k  +  *k+l(i))  1  (k+1)  (2.27) 


where  x^d)  is  the  value  of  the  i-th  point  (t  =  0.12i  sec)  of  the  (k+l)-th 

run. 


2.  The  ensemble  mean-square  x.+^(i)  was  computed  similarly  from 
Xk+l(1)  =  (\0)'k  +  x^+1(i))  /  (k+1) 


(2.28) 


and  the  unbiased  estimate  of  the  ensemble  variance,  V^d),  was  then  found  in 
the  usual  manner 

Vk(1)  =  (\(i)  -  \(i))'  k/(k-D  (2.29) 

This,  of  course,  resulted  in  the  ensemble  mean  and  ensemble  standard  deviation, 
viz. , 

N 

x(t)  =  £  ^  xk(t)  (2.30) 

k=i 


N 

Nd  S  [Xk(t)  "  ^(t)]2 

k=l 

The  resulting  time  histories  are  plotted  as  follows: 

•  Figures  11  -  12  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 


1/2 


(2.31) 


ax(t)  = 


Static-G,  0-peak. 


4  1 


•  Figures  13-14  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 

Static-G,  1-peak 

•  K inures  15  -  lb  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 

Static-G,  5- peak. 

•  Figures  17-18  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 

Dynamic-G,  0-peak. 

•  Figures  19  -  2u  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 

Dynamic-G,  1-peak. 

•  Figures  21-22  :  Mean  and  SD  of  Longitudinal  Tracking  Error, 

Dynamic-G,  5-peak. 

In  Figures  23  -  29  only  the  initial  60-second  period  of  the  1-peak 
profile  is  considered. 

•  Figures  23-24  :  Mean  and  SD  of  Lateral  Tracking  Error,  Static-G. 

•  Figures  25-26  :  Mean  and  SD  of  Lateral  Tracking  Error,  Dynamic-G. 

•  Figure  27  :  Commanded  G  Mean,  Static-G. 

z 

•  Figure  28  :  Commanded  G  Mean,  Dynamic-G. 

z 

•  Figure  29  :  Subject-Attained  G  Mean,  Dynamic-G., 

Notice  that  not  all  variables  are  plotted  here.  This  is  done  mainly  for  practical 

reasons.  For  example,  the  Commanded  G  for  the  5-peak  trajectory  was  an  exact  replica 

z 

of  the  1-peak  case.  The  SD  of  this  variable  is  of  no  great  consequence  since  it  is 
virtually  zero.  Also,  the  stick  inputs  in  both  axes  were  unreadable  from  the  supplied 
tapes.  These  variables,  however,  are  of  minor  interest  in  the  modeling  process. 

The  next  step  in  the  processing  of  the  ensemble  data  was  statistical  analysis. 

Since  the  air-to-air  experiments  were  performed  under  two  experimental  conditions 
(Static-G  vs.  Dynamic-G),  statistical  tests  had  to  be  employed  in  order  to  detect  any 


pig.  12  STANDARD  DEVIATION,  0  PEAK,  STATIC  G 
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Fiq.  23  MEAN  HEADING  TRACKING  ERROR,  1  PEAK,  STATIC  G 


Fig.  24  STANDARD  DEVIATION  HEADING  ERROR,  1  PEAK,  STATIC  G 


►9 


Fiq.  25  MEAN  HEADING  TRACKING  ERROR,  1  PEAK,  DYNAMIC  G 


Fiq.  26  STANDARD  DEVIATION  HEADING  ERROR,  1  PEAK,  DYNAMIC  G 
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F'<|.  <?9  SUBJECT  ATTAINED  G  MEAN,  DYNAMIC  G 
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statistical  difference  among  the  ensemble  means  and  variances.  For  convenience,  only 
the  first  60  seconds  of  the  1-peak  trajectory  are  considered,  but  the  results  can  be 
generalized  for  the  other  target  profiles.  The  following  hypotheses  were  examined: 

(11  I.ongitudinal  tracking  error  ensemble  means  are  equal  under  both  experimental 
conditions.  (Welch-test  )  (Figure  3U1 . 

(2)  LongitudinaL  tracking  error  variances  are  equal  under  both  experimental 
conditions  (F-test)  (Figure  31). 

(3)  Same  as  (1)  for  the  lateral  dalu  (Figure  32). 

(4)  Same  as  (2)  for  the  lateral  data  (Figure  33). 

Application  of  the  statistical  tests  and  inspection  of  Figures  30  -  33  unveils 
the  following: 

a.  We  reject  hypothesis  (1).  A  comparison  of  Figure  30  with  Figures  7  and 

29  (target's  G-profile  (GT)  and  attacker’s  attained  G  (G  ),  respectively) 

shows  that  there  is  a  significant  difference  (P<.05)  between  the  means  during 

active  G  (12-18th  second  of  the  run)  and  during  high  G  (around  the  30th 
L  A 

second  of  the  run)  periods. 

b.  We  reject  hypothesis  (2).  It  is  evident  that  the  variances  are  significantly 
different  (P<.05)  during  a  part  of  the  active  GT  period  (ll-16th  second)  and 
during  most  of  the  period  that  immediately  follows  the  peak  G  interval. 

c.  We  fail  to  reject  hypothesis  (3).  This  result  is  obvious  since  there  is  no 
target  input  in  the  horizontal  axis.  Also,  the  acceleration  forces  are 
perpendicular  to  this  plane. 

•j* 

In  order  to  test  the  null  hypothesis  of  equality  of  means,  one  would  normally  use 
the  t-test  only  when  equality  of  variances  is  assumed.  In  our  case,  however,  the 
null  hypothesis  does  not  assume  to  deal  with  the  same  population,  since  the  variances 
prove  to  be  significantly  different  (see  (2)  and  (4)).  In  such  a  case  Welch  test  is 
usually  employed.  In  the  event,  both  t  and  Welch  tests  gave  almost  identical  results. 
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il.  We  reject  hypothesis  (A)  (P<.05).  A  significant  difference  between  the 
lateral  tracking  error  variances  is  evident  (Figure  33). 

An  inspection  of  the  aggregated  data  indicates  the  following. 

1.  In  general,  the  quality  of  the  data  is  excellent.  There  is  a  sufficient 

number  of  runs  to  form  a  "good"  ensemble  average.  Moreover,  unlike  the 
previous  High  Acceleration  Cockpit  (HAC)  studies,  the  averages  are  relatively 
smooth,  tend  to  he  more  believable  statistically,  and  provide  a  good  corner¬ 
stone  for  C~modeling  efforts. 

.  The  G  vs.  no-G  comparisons  of  the  longitudinal  tracking  errors  and  of  the 

standard  deviations  (e^(t)  and  o  i  t ),  respectively)  show  that  the  magnitudes 

of  e  (t)  are  not  very  different.  There  is,  however,  an  increase  in  0  (t) 

p  ep 

for  the  G-present  versus  the  static  runs.  This  is  in  agreement  with  past 
acceleration  research  wherein  performance  decrements  are  observed  with  sustained 


G-stress . 

3.  It  is  noted  that  the  effective  time-constant  of  e  (t)  in  going  from  the  last 

P 

7  G  peak  to  the  A  G  tracking  phase  is  smaller  in  the  static  vs.  the  dynamic 

runs.  Thus,  the  time  required  for  e^(t)  ^  0  is  longer  with  the  G-stress 

present.  Similar  time-constant  differences  can  be  observed  in  the  0  (t). 

ep 

I'h  i  s  slower  pilot  response  is  again  compatible  with  previous  G-stress  results. 

A.  The  static  tracking  e  data  for  the  1-  or  5-peak  trajectory  show  a  tendency 

P 

for  the  error  to  return  towards  zero  immediately  following  a  target  pull-up 
maneuver.  Thus,  the  peak  (negative)  error  exhibits  a  positive  slope  tendency. 
5.  The  same  e^(t)  tracking  data  for  the  G-stress  case  show  a  diametrically 
opposed  trend.'  During  the  peak  error  period  there  is  a  further  lagging 
it*  nd  i:i  the  subject  's  response,  i.e.  the  error  decreases.  This  tendency  was 
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noted  in  the  earlier  HAC  data,  but  could  not  be  confirmed  as  being  statisti¬ 


cally  significant.  Indeed,  it  was  suggested  by  some  that  this  tendency  was 
due  to  pilot  anticipation  of  the  target’s  upcoming  pitch-rate  reversal. 

This  explanation  is  not  borne  cut  by  the  new  data,  however.  As  we  follow  the 
ensemble  average  of  e^(t)  through  the  5  peaks  this  lagging  tendency  becomes 
diminished.  If  there  were  pilot  anticipation  it  should  have  shown  up  through¬ 
out  the  run,  nossibly  becoming  more  pronounced  on  subsequent  peaks.  Instead 
we  see  the  opposite  trend.  It  appears  that  the  person  may  becoming  "acclimated" 
to  the  repetition  of  G-stress  and  local  performance  is  improving.  Put  another 
way,  the  degrading  effects  of  the  G-stress  arc  possible  being  absorbed — the 
pilot's  G  threshold  might  actually  increase. 

6.  The  results  do  not  appear  to  be  tainted  by  the  lateral  axis  control  tasks. 

On  the  contrary:  a  valuable  insight  into  the  effects  of  G-stress  on  pilot 

performance  is  gained  by  the  inclusion  of  this  secondary  task.  The  lateral 

tracking  error  means  (e^Ct))  are  not  substantial,  which  is  natural  as  there 

is  no  target  input  in  this  axis.  The  standard  deviation  (a  (t))  magnitudes, 

however,  are  worth  noticing.  It  is  readily  evident  that  O  (t)  in  the  G- 

er 

stress  condition  is  much  larger.  This  is  a  testimony  to  the  thesis  that 
there  is  a  substantial  performance  degradation  when  subjected  to  sustained 
G-stress . 

7.  It  appears  that  the  Attained  G  profile  is  a  delayed  and  filtered  version  of 
the  Commanded  G.  This  is  because  of  the  slow  dynamics  of  the  l)ES,  which 

can  be  approximated  by  a  first  order  low-pass  filter,  with  time  constant  ~  4 . 5 
sec . 
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2 .  3  Modeling  the  Longitudinal  Tracking  Task 
2.3_.  1  Appli cation__of^  the  Optimal  Control  Model 

The  Optimal  Control  Model  (OCM) ,  modified  to  treat  deterministic  target  motion, 
assumes  the  system  dynamics 


x  (t)  =  A  x  (t)  +  b  6  (t)  +  F  z(t) 
-o  o-o  o  A  — o 


y(t)  =  C  x  (t) 
0-0 


(2.32) 


where  A  ,b  ,F  and  C  are  the  pertinent  system  matrices,  x(t)  is  the  system  state 
vector,  6  (t)  is  the  elevator  deflection  and  y_( t )  are  the  observations  available  to 
the  human  operator  (HO)  . 

One  of  the  difficult  aspects  of  the  modeling  is  the  choice  of  state  variables. 
The  criteria  are 

1.  a  and  q  for  both  aircraft  should  be  retained  as  states. 

2.  The  equations  should  be  simple. 

3.  Other  variables  such  as  G  ,  e  ,  r,  etc.,  should  either  be  states  or  linear 

z  p  ’ 

combinations  of  states. 

4.  No  absolute  (i.e.  inertial)  angles  may  be  states 


5. 


The  rates  e  and  r  should  be  combinations  of  states. 
P 


The  last  criterion  indicates  that  a  relative  coordinate  system  is  to  be  used. 

After  considering  numerous  possibilities,  the  following  choices  were  made 


This  condition  is  needed  for  subsequent  pilot  modeling. 
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V 


xi  =  qT 


x2  =  aT 


X3  =  qA 


X4  =  aA 


x.  =  0A  -  0 
5  A  T 


x,  =  e  =  0.  - 
6  p  A  T 


(2.33) 


Target  dynamics 


Attacker  dynamics 


(2.34) 


Differentiating  these  quantities  gives 


•  T  T  T 

x  =  M  x-  +  M  X,  +  Mj  6_ 
1  a  2  q  1  6  T 


x2  =  Za  x2  +  X1 

x,  =  x.  +  x-  +  6 

3  a  4  q  3  6  A 


x  ,  =  Z  x,  +  x- 
4  a  4  3 


X5  x3  '  X1 


x6  =  x3  +  D  (X2  “  x4  +  x5> 


The  pilot  observes  the  error  e  (t)  and  aspect  angle  r(t),  plus  their  rates  of  change. 

P 

Therefore 


y.  =  e  =  x, 

1  p  6 


y2  =  ep  =  X3  +  D  (X2  '  x4  +  x5> 


y.  *  r  ■  9_  -  JL  ■  x,  -  x- 
7  3  T  t  6  5 


(2.35) 


y4  “  *  “  x6  '  X5  “  X1  +  5  ’  (X2  '  X4  +  V 
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Also  note  that 


G 

z 


V  ^  ^ 
a  A 


(2.36) 


In  order  to  simplify  the  modeling  process  we  sought  a  reduced-order  version  of 

the  model.  Therefore  we  assumed  that  the  attacker  and  the  target  aircraft  flight 

path  angles,  y^  and  y  ,  equal  0A  and  0  ,  respectively.  Note  that  in  this  situation 

the  state  variables  x.  and  x,  can  no  longer  be  regarded  as  or  and  a_,  since,  strictly 

2  4  A  L 

speaking,  or  =  i  =  0.  However,  tin-  2::d  order  short  period  attacker  dynamics  remain 

a  r 

as  in  Equation  (2.34).  The  assumption  y  =  0  modifies  only  the  equation  for  x . , 

A  A  D 


X6  =  0A  -  :'t  =  X3  "  D  (0T  '  V 


(2.37) 


X3  +  D  X5 


The  result  of  the  target  simplification  is  that  state  x^  does  not  enter  the  system 
equations,  except  for  the  target  dynamics.  If  the  explicit  2nd  -  order  response  of 
the  target  is  not  important  for  modeling,  then  it  is  convenient  to  replace  the  equation 
for  Xj  with 

x^  =  z(t)  (2 . 38) 

where  x^  =  qT(=yT)  and  z(t)  is  the  target's  commanded  flight  path  angle  acceleration 
(i.e.  pitching  moment). 

Taking  these  simplifications  into  account  and  redefining  states  to  account  for 
omission  of  >  we  find 
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These  modeling  assumption  are  in  accordance  with  the  simulation  setup  at  AMRL. 

There  is  a  subtlety  involving  the  vertical  acceleration  C^:  while  it  is  always 

true  that  G  scales  with  it  will  not  scale  with  a,  in  this  case.  But  since 

z  A  A 

Y.  =  q .  ,  the  vertical  G  force  should  be 
A  A  z 


G 

z 


(2.40) 


Therefore,  the  vertical  acceleration  of  the  target  aircraft  — in  g  units —  is 


v  •  V  V 

GT(t)  -  g  Yt  =  -  qT  =  -  Xj. 


(2.41) 


and  the  vertical  acceleration  commanded  by  the  attacker  is  given  by 


GA(t)  "  I  V‘>  +  1  =  g  X2(t)  +  1 


(2.42) 


Because  of  distributed  dynamic  lags,  nonlinear  effects,  etc.,  the  actual  accelera¬ 
tion,  G  ,  attained  on  the  centrifuge,  and  hence  that  felt  by  the  pilot  is  approxi¬ 
mately 


6  I 


G.  =  C(s)  G^,  (2.43) 

A  A 

where  C(s)  =  -  approximates  the  DES '  transfer  function. 

1+4 . os 

2.3.2  G-Submodel  for  Longitudinal  Tracking  Task 

The  OCM  assumes  that  the  well-trained  and  well-motivated  human  operator  adopts 
an  optimal  control  strategy,  subject  to  his  various  psycho-physiological  limitations. 

The  task  of  correctly  specifying  these  limitations  along  with  the  resultant  control 
strategy  is  not  always  easy  since  they  depend  heavily  on  environmental  and  experimental 
conditions.  In  the  case  at  hand,  the  environmental  condition  is  the  level  of  G-stress 
the  human  is  subjected  to.  Below  we  describe  the  major  human  limitations  that  are  of 
concern  in  our  modeling  efforts. 

Operator  Time  Delay.  The  various  internal  time  delays  associated  with  visual  and 
central  processing  are  represented  by  an  equivalent  time  delay,  T.  The  sensitivity 
analysis  that  preceded  the  modeling  phase  had  not  uncovered  any  particular  trends 
in  the  model  predictions.  It  proved  to  be  largly  insensitive  to  this  parameter.  There¬ 
fore,  the  nominal  value  of  T  =  0.25  seconds  was  chosen  for  the  Static-as  well  as  for 
the  Dynamic-G  condition. 

Threshold  Effects:  Nonlinear  threshold  effects  associated  with  human  information 
processing  are  the  visual  and  the  indifference  thresholds  (not  distinguished  by  the 
model).  We  therefore  associate  with  each  displayed  variable  y^(t)  a  visual  indifference 
threshold  level,  Thus,  the  (delayed  and  noisy)  perceived  signal  is 

ypl(t)  =  Fityi(t-T)]  +  vyl(t-T)  i“l , 2 , 3,4  (2.44) 

where  v  .(t)  is  assumed  to  be  a  white-Gaussian  observation  noise  (see  the  following 

yi 

paragraph),  and  the  threshold  function  F^(>)  is 
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2* 


0  -a.  '  V.  '  n t  1  =  1, ...  ,4 

yi  "  al  3i  -  y i  (2.A5) 

Since  v(t)  is  assumed  Gaussian,  we  apply  here  the  concept  of  the  Random- Input- 
Descr  ib  ing-Funct  ion  (RIDF)  to  statistically  linearize  F(-)  [43].  The  approximated 

perceived  observation  is  given  therefore  by 

v  .(t)  ~  N  v.(t-T)  +  v  .(t-T)  i= 1,2, 3,4  (2.46) 

pi  i'  l  yi 

where  N  =  N.(a.,v.,0  .)  is  the  i-th  effective  display  gain,  y.  being  E(y  )  and 
i  l  i  ■  i  vi  l  i 

0  .  being  v  's  standard  deviation.  The  nominal  values  of  a  ,  a  ,  a  were  determined 
vi  i  234 

from  the  experimental  setup  and  the  displayed  target  size; 

a^  =  2mrad/sec,  a^  =  50  mrad  a^  =  25  mrad/sec 

In  both  experimental  conditions  these  values  remained  unchanged.  The  choice  of  the 
indifference  threshold  on  the  tracking  error,  a^,  was  a  more  subtle  task,  and  it  is 
discussed  later  in  the  section.  It  will  be  seen  that  a^  is  heavily  dependent  on 
tfie  (,-stress  level  the  tracker  is  subjected  to. 

Observation  and  Motor  Noise:  The  various  sources  of  inherent  human  randomness  are 
represented  by  observation  noise,  v  ,  and  motor  noise,  v^.  These  zero-mean  white 
noises  mainly  account  for  errors  in  observing  displayed  variables  and  in  generating 
control  signals.  The  observation  and  motor  noise  covariances  are 
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i=l . A 


(2. 47) 


V  . 
yi 


(t) 


Vg(t)  -  Tip^  •  cov[6A(t)] 

where 


(2.48) 


p  (pr)  =  observation  (motor)  noise  ratio  coefficient 

y  o 

f.  =  fractional  attention  allocation  to  indicator  i 

l 

The  nominal  values  for  the  noise/signal  coefficients  were 

pyi  =  -21.5  dB  i=l,2, 3,4 

P6  =  -18.0  dB 

These  values  were  applied  to  the  Static-G  condition.  In  the  Dynamic-G,  p  remained 
unchanged  but  the  motor-noise  ratio  coefficient,  pg,  had  to  be  increased  by  ~  40% 

(to  -16.5  dB)  to  match  the  empirical  data,  yielding  a  40%  increase  in  the  motor 
noise,  V^(t). 

Since  there  are  four  observation  channels  —  two  for  the  longitudinal  axis 
(e^(t)  and  r(t)),  and  two  for  the  lateral  (ef(t)  and  ( t ) )  —  the  fractional  atten¬ 
tions  are  assumed  to  be  evenly  divided  among  them.  The  result  is  f  =  1/4,  which 
increases  the  effective  observation  noise  by  6  dB. 

One  may  notice  that  in  all  the  covariance  propagation  equations  of  Section  1.4, 
which  involve  =  diag  [V^],  (Equations  (1.18,  (1.30)  -  (1.32)),  only  C  C 
appears.  If  we  let 

N  =  diag  (N^)  (2.49) 

and 

v  >'r  i  =  col  [v  .  (t)l,  (2.50) 

v  vi 
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then  Equation  (2.46)  becomes 


y  (t)  =  NCx(t-T)  +  v  (t)  (2.51) 

p  -  -y 

It  can  be  seen,  therefore,  that  the  matrix  C  has  to  be  modified  to  NC.  We  preferred, 
however,  to  modify  viz., 

V_1  -*■  N'V_1N  .  (2.52) 

y  y 

The  mathematical  equivalent  is, 


ypi(t)  =VC~T)  + 


\i(t) 


(2.53) 


-Pyi 


lyj(t)  +  a2n) 


(2.54) 


i  i 

which  corresponds  to  Equations  (1.4)— (1.5).  This  interpretation  of  the  visual  channel 


is  very  significant.  The  net  effect  of  the  nonlinearity-equivalent  gain  is  to 
increase  V^^t),  since  from  RIDF  theory,  0  <  <  1.  In  fact,  decreases  with 

increasing  threshold  level  a^.  It  will  be  seen  that  the  tracking  error  indifference 
threshold,  a^ ,  increases  with  G,  thus  decreasing  N^(a^),  which  in  turn  increases  the 
observation  noise  covariance  of  the  tracking  error,  Vg^(t)  =  V^(t). 


Neur  motor  Time  Constant:  Although  it  has  been  found  that  the  model  is  quite 

insensitive  to  variations  in  the  neuromotor  time  constant  t„,  we  have  chosen  to 

N 

adopt  the  values  of  »  .18  sec  for  the  Static-  and  ■  .2  sec  for  the  Dynamic-G 
condition.  This  parameter,  which  governs  the  first-order  lag  in  the  human  control 
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mechanism,  has  proven  to  be  a  key  parameter  in  the  steady-state  G-model  (Chapter  III). 
A  similar  difference  in  the  time-constant  values  between  the  two  experimental  condi¬ 
tion  was  identified  there. 


Control  Cains:  In  computing  the  optimal  gains  adopted  by  the  human,  we  minimize  the 
following  cost  functional, 


J (6a)  =  £im  - 
A  T  T 

i-*oo 


[qe  ej(t)  +  ^(t)  +  g  6A(t)  ]  dt 


(2.55) 


where  q^,  q^  and  g  are  the  (constant)  weighting  coefficients.  The  control  rate 
weighting  g  is  completely  determined  by  the  choice  of  q  ,  but  q^  and  q^  are  subjective 
and  represent  the  operators  control  strategy  and  error  tradeoffs.  We  assume  q^=q^=.l 
in  both  conditions  since  in  the  longitudinal  tracking  task  the  subjective  importance 
of  the  pitch  error  rate  is  no  lesser  than  that  of  the  pitch  error. 

In  the  Repetitive-G  experiments,  there  were  two  major  factors  that  influenced 
the  performance  of  the  subjects: 

(1)  The  specific  instructions  and  training  that  guided  the  subjects  in  their 
task; 

(2)  The  attained  positive  acceleration  (G^)  which  resulted  from  the  commanded 

c 

input  signal  (G^) . 

It  is  clear  that  coaching  the  subjects  to  adopt  only  certain  prespecified  control 
actions,  and/or  constraining  their  permissible  command  levels,  should  profoundly 
affect  tracking  performance.  This  was  demonstrated  by  Brown  and  Collins  [26]  in 
experiments  whereby  subjects  were  required  to  perform  a  tracking  task  similar  to 
ours.  Prior  to  certain  runs  the  subjects  were  instructed  to  fly  the  simulated  air- 
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craft  in  a  coordinated  fashion,  i.e.,  to  maintain  lateral  acceleration  at  a  minimum. 
The  results  proved  that  there  was  a  significant  difference  in  performance  between 
these  runs  and  those  not  preceded  by  the  specific  instructions. 

The  effect  of  such  instructions  is  reflected  in  the  human's  control  strategy, 
i.e.  in  the  weightings  in  the  cost  functional  J(-)*  In  our  experiments,  the  subjects 
were  instructed  (and  trained)  to  track  well,  i.e.,  to  minimize  the  tracking  error, 
hut  only  during  the  last  stage  of  the  run,  namely,  from  the  onset  of  the  4  0  stress 
period.  The  subjects  were  scored  only  during  the  final  90  second  period  of  the  4  C, 
plateau.  Preceding  the  4  G  period,  the  subjects  were  instructed  to  merely  "keep  the 
target  on  the  screen." 

One  might  well  conclude  from  this  experimental  setup  that  the  subjects'  error 
indifference  threshold  level,  a^,  was  changing  during  the  tracking  period,  and  that  it 
had  a  profound  effect  on  the  relative  error /error-rate  weightings  in  the  computation 
of  the  optimal  gains.  One  plausible  error  threshold  time  history  is  shown  in  Figure 
34.  Figure  34  refers  only  to  the  1-peak  trajectory  but  similar  threshold  time- 

histories  mav  be  used  for  the  other  trajectories.  During  the  first  tracking  interval, 
a^(t)  =  thj  =  maximum  angular  error  that  corresponds  to  the  graphics'  total  screen 
size  =  150  mrad .  During  the  G^,  =*  4  G  plateau  a  (t)  =  thQ  =  minimal  threshold  level 
inherent  to  the  human,  and  dependent  on  the  size  of  the  pipper  image  simulated  on  the 
screen  (=  20  mrad).  We  assume  that  the  threshold  changes  gradually  around  t  =  30 
sec,  where  G^.  decreases  from  7  to  4  G.  These  assumption  were  applied,  regardless  of 
ttie  experimental  condition  (i.e.,  for  both  Static-  and  Dynamic-G) . 
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With  Che  above  In  mind,  coupled  with  Equations  (2.46)  and  (2.55),  we  replace  the 
displayed  pitch  error  ep(t)  in  the  cost  functional  J(*)  with  the  perceived  error, 
namely. 


where 


*P<t)  -  "i !•!<*>.  %<*>•  °.p(t)J  «p(t) 


ep(t)  -  y1(t) 


°ep(t)  *  °yl(t) 


Thus,  the  cost  functional  to  be  minimized  is  now 

T. 


(2.56) 


(2.57) 


J,<V 


’  T~  T  E{ J  lqeNl  ep(t)  +  qe  «p  +  8  *A>  dtj 

"O 


(2.58) 
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The  error  rate  display  gain  is  not  included  in  J'(-)  since  we  assume  that  the  error 
rate  threshold  level,  a ^  is  negligible  (=  2  mrad/sec)  and  time-invariant.  This 
approach  of  weightings  modification  in  the  cost  functional  has  been  originally 
suggested  by  Kleinman  and  Baron  [A3],  and  it  extends  the  preliminary  modeling  efforts 
of  Korn  et  al  [AA]. 

One  may  argue  that  this  approach  has  an  inherent  serious  drawback.  Since  a^  is 
now  time  dependent,  the  error  weighting  would  vary  in  time  which  in  turn  implies  that 
the  optimal  control  gains  would  have  to  be  computed  continually.  This  would  not  be 
feasible.  However,  it  has  been  found  that  only  the  error  (x,.)  gain  depends  strongly 
on  the  error  weighting,  and  the  remaining  gains  are  virtually  unaffected.  With  this 
fact,  a  functional  relationship  between  the  error  gain  and  the  error  weighting  has 
been  found,  and  successfully  applied. 

The  indifference  threshold  on  the  observed  tracking  error  is  further  affected  by 
the  G-stress.  We  postulate  that  the  human  would  tend  to  decrease  the  importance  he 
assigns  to  the  tracking  performance  when  subjected  to  sustained  physical  stress. 

From  Figures  17  -  22  it  is  apparent  that  the  subjects  exhibited  a  lagging  tendency, 
which  increased  during  the  peak-G  periods.  The  hypothesis  is,  therefore,  that  the 
acceleration  stress  would  tend  to  increase  the  tracker's  indifference  threshold 
associated  with  the  tracking  error.  Alternately,  he  would  become  less  indifferent 
to  his  own  psycho-physiological  response. 

To  incorporate  this  threshold  effect  in  our  G-model  we  adopted  an  ad-hoc  formula, 
such  that , 

al(t’  CA>  =  al^C)  l1  +  pg-max2  (0»  Ga  “  ^mln^  ^  (2.59) 
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where  G  =  G  (t)  is  the  subject  attained  G  (Equation  (2.43))  and  G  .  =  4g  =  minimal 

A  A  z  min 

(■  level,  below  which  human  performance  is  assumed  not  to  be  affected  significantly. 

(We  assume,  et  course,  that  the  subject  is  wearing  a  fl-suit  j  .  The  value  of  the  parameter 


that  gave  the  best  data-model  fits  was  p  =  .5. 

8  8 

The  effect  of  the  increase  in  a,  under  G  -stress  is  twofold: 

1  z 

1.  The  equivalent  display  gain,  N^,  decreases  under  G-stress,  further  reducing 
the  effective  weighting  on  the  tracking  error,  while  the  error-rate  weighting 


remains  unchanged.  This,  in  turn,  reduces  the  control  gain  on  the  error  that 


is  computed  from  Equation  (2.58). 


2.  The  decrease  of  increases  the  effective  observation  noise  covariance 

associated  with  the  pitch  error  (V  (t))  as  can  be  observed  in  Equation 

ep 

(2.54). 


2.3.3  Modeling  Results  of  the  Longitudinal  Tracking  Task 


A  good  performance  model  is  one  that  satisfactorily  validates  experimental  data 
and  can  be  applied  reliably  in  a  predictive  mode.  Moreover,  it  should  be  applicable 
to  a  wide  range  of  manual  control  tasks.  The  OCM  has  proven  its  ability  to  model 
the  air-to-air  compensatory  tracking  task  as  demonstrated  by  the  excellent  data 
matches  that  follow. 

In  order  to  demonstrate  the  model  predictions  of  the  pertinent  variables,  we 
consider  the  1-peak  trajectory  as  a  representative  case.  The  0-peak  profile  is  of 
no  real  interest  because  of  a  limited  target  maneuver  and  the  low  G  levels,  and  the 
5-peak  profile  is  basically  a  replica  of  the  1-peak  as  can  be  observed  from  the 
ensemble  data.  Moreover,  most  of  the  final  period  of  constant  4  G  can  be  discarded 
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since  it  is  a  period  of  zero  q  (z(t)  =  0).  Therefore,  it  is  sufficient  to  consider 
the  first  60-second  period  of  the  1-peak  profile,  in  order  to  illustrate  our  modeling 
resul ts . 

In  the  following  discussion  we  present  the  model  prediction  time  histories. 

These  are  followed  by  figures  which  show  the  model-data  comparisons. 


• 

F  gu  rt  . 

i  35  - 

36: 

Mean  and  SD  of 

tracking  error, 

Static-G. 

• 

Figures 

i  37  - 

38: 

Mean  and  SD  of 

tracking  error. 

Dynamic  G 

« 

Figure 

39: 

Commanded-G,  Dynamic-G 

condition. 

• 

Figure 

4u : 

Attained-G  profile. 

Comparisons : 

•  figure  4  1  :  fracking  error  mean,  Static-G. 

•  Figure  4 1  •  fracking  error  SD,  Static-G. 

•  Figure  43:  fracking  error  mean,  Dynamic-C. 

•  Figure  *4:  Tracking  error  SD,  Dvnamic-G. 

•  Figure  -t5:  Commanded  G,  Static-G  condition. 

•  Figure  46:  Commanded  G,  Dynamic-G  condition. 

2.3.4  Discuss  ion 

The  data-model  fits  are  excellent.  It  should  be  emphasized  that  the  base-line 
OCM  parameters  were  unperturbed  from  their  nominal  values  in  the  no-stress  condition. 
Only  in  the  G-stress  situation  was  it  necessary  to  formulate  a  G-submodel.  To 
recapitulate  the  G-stress  performance  model,  the  model  predictions  have  to  be 
expounded : 
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Fig.  41  MODEL-VS-DATA,  PITCH  ERROR  MEAN,  STATIC  G 


pPJW 
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CH  ERROR  MEAN,  DYNAMIC  G 


Fig.  44  MODEL- VS- DATA,  PITCH  ERROR  SD,  DYNAMIC  G 


MODEL- Vr> -DATA,  COMMANDED  G  ,  STATIC  G  CONDITION 


8> 


Fig.  46  MODEL- VS- DATA,  COMMANDED  G  ,  DYNAMIC  G  CONDITION 


(1)  One  can  observe  a  slower  recovery  of  the  tracking  error  during  the  transition 
period  between  the  7 G peak  to  the  4 G final  plateau,  when  under  G-stress. 

This  model  prediction  is  attributed  to  the  increase  in  the  indifference 
threshold  on  the  tracking  error,  a^,  in  the  Dynamic-G  condition.  The  major 
effect  of  increasing  a^  is  to  decrease  N^a^),  which  in  turn  decreases  the 
relative  weighting  on  the  tracking  error  in  Equation  (2.58).  The  resulting 
feedback  gain  on  the  tracking  error  is  lower  in  magnitude,  allowing  the 
introduction  of  larger  tracking  errors,  thus  permitting  a  slow  error  recovery 
Notice  (figure  40)  that  the  peak  Attained-G  level  occurs  during  this  transi 
tion  interval. 

(2)  A  secondary  effect  of  the  increase  in  a^  under  G-stress,  is  to  effect  higher 
observation  noise  V  (t)  (Equation  (2.54)). 

(1)  Another  indication  of  performance  degradation  under  G-stress  is  larger 

standard  deviations  of  tracking  error.  This  is  demonstrated  in  figure  38. 

The  model  predicts  a  larger  and  persisting  tracking  error  variance  in  the 
Dynamic-G  condition,  which  is  attributable  to  higher  motor-noise.  Indeed, 
a  40%  higher  motor  noise  (p^  =  -16.5  dB  in  Dynamic-versus  =  -18  dB  in 
Static-G)  was  necessary  to  fit  the  ensemble  tracking  error  SD.  (From  the 
covariance  propagation  equations  of  Section  1.4,  one  can  see  that  increasing 
the  motor-and/or  observation  noise  results  in  higher  variances.) 

The  modeling  assumptions  prove  to  be  in  excellent  agreement  with  past  G-stress 
performance  research,  as  summarized  in  Section  1.2.  Virtually  all  of  these  efforts 
conclude  that  the  performance  decrements  are  ascribed  to  larger  motor  and  observation 
noises,  as  well  as  to  such  psychomotor  response  as  tolerating  larger  tracking  errors 
when  operating  under  sustained  acceleration  stress. 


2.4.1  Application  of  the  Optimal  Control  Model 


Similar  to  the  longitudinal  tracking  dynamics,  the  lateral  task  has  to  be 
expressed  in  the  state  space. 


x  (t)  =  A  x  (t)  +  b  6  (t) 
— o  o-o  — o  a 


y(t}  =  C^t) 


(2.60) 


where  A  ,  b  and  C  are  the  system  matrices  of  the  lateral  dynamics,  x  (t)  is  the 
o—oo  — o 

state  and  <5^(0  is  the  aileron  deflection.  Notice  the  absence  of  the  driving  terms 
w  (t)  and  z(t).  In  accordance  with  the  simulated  roll  dynamics  as  represented  by 
Gr(s)  in  Section  2.1.5,  we  define  the  state  as 


x^  =  e^(t)  =  lateral  tracking  error 


X2  =  er(t)  =  D  * 


x3  =  K(t)  <l> 


x4  =  K(t)  i|/ 


(2.61) 


It  should  be  noted  that  these  definitions  are  valid  only  when  K(t)  is  constant  or 
slowly  varying.  This  has  been  the  usual  assumption  when  considering  time-varying 
system  parameters  in  pilot  modeling  efforts. 


This  4-th  order  system  can  be  reduced  to  third  order  by  combining  the  first 

order  Lag  of  Equation  (2.26),  (t.  s+1)  \  .  •  *>  the  first  order  lag  (t.,  s+1)  *  which 

q)  N 

is  introduced  by  the  neuromotor  dynamics  in  the  OCM.  Thus  we  set 
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T  =  T  +  T  -  T  + 

N  ,  equiv  N  <f>  N  5.55 


in  the  subsequent  modeling  process,  and  the  system  then  becomes 


X1  X2 


x^  =  x^ 


x.  =  K(t)6  (t) 
3  a 


(2.62) 


(2.63) 


where  i  (t)  = 


obta  in , 


ip,  and  x^,  x3  are  as  defined  previously.  In  matrix  form,  we 


0  10 


x  =  OOlx+O  K(t)  6  (t) 
— o  ~ o  a 


(2.64) 


0  0  0 


Note  that  K(t)  enters  only  as  a  scalar  multiplier  on  the  control  6^(0.  This 

means  that  the  control  6'(t)  can  be  computed  for  K(t)  =  1  and  then  adjusted  for 

a 

arbitrary  K(t)  via  inverse  scaling.  It  will  be  seen  that  K(t)  has  an  effect  on 
the  pilot's  additive  motor  noise.  For  notational  purposed  we  define  K(t)6a(t)=u(t) 
The  observations  that  are  available  to  the  pilot  to  minimize  lateral  tracking 
error  are 


yl  *  X1 


y2  =  X2 


=  tracking  error  e(t) 


=  error  rate,  e(t) 


(2.65) 


y 3  =  x3/K(t)  = 


'  “A t) 


relative  roll  angle,  <p 


roll  rate. 


H3 


The  observation  of  roll  angle  is  critical  for  the  control  task.  If  this  inform- 
ation  is  absent  it  is  virtually  impossible  for  a  subject  to  control  the  K/s  system. 

The  rule-of-thumb  control  logic  of  "keep  your  wings  aligned  with  those  of  the 
target  aircraft"  is  testimony  to  this  fact.  The  information  y^  is  not  very  impor¬ 
tant  here  inasmuch  as  the  control  signal  u(t)  is  essentially  "known."  It  would  be 
of  importance  if  the  target  aircraft  were  free  to  move  laterally. 

2.4.2  G-Submodel  for  Lateral  Tracking  Task 

Similar  to  the  longitudinal  tracking  task,  our  modeling  approach  was  first  to 
isolate  the  G-stress-dependent  OCM  parameters  and  cost  functional  components,  and 
then  to  determine  the  structure  of  that  dependency.  Below  we  list  all  OCM  parameters 
and  indicate  how  they  are  affected  by  acceleration  stress: 

Threshold  Effects:  The  threshold  nonlinearities  were  applied  here  in  a  fashion 
similar  to  that  of  Section  2.3.2.  The  nominal  indifference  threshold  values  on 
y^,  y3  and  y^  were  selected  from  simulation  considerations  as  =  2  mrad/sec, 
a^  *  5  rarad,  3  2.5  mrad/sec.  These  values  remained  unchanged  in  the  dynamic  condi¬ 
tion.  The  threshold  level  on  the  observed  lateral  tracking  error  was  chosen  here,  as 
in  Section  2.3.2,  to  reflect  the  pilot-adopted  tracking  strategy.  In  the  nominal 
static  condition  we  select 

t<tg  a^(t)  =  thj  3  35  mrad  »  maximum  angular  error  allowable 
t>tg  a^(t)  =  thQ  3  25  mrad  =  pipper  radius. 

In  the  G-stress  condition  th.  *  50  mrad,  and  the  value  of  th  is  unchanged.  In  the 

1  o 

presence  of  G-stress,  there  is  assumed  to  be  a  further  increase  in  the  model's 


\ 
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indifference  threshold,  viz; 


al(t,<V  =  al(t)  ^  +  pg' 1113x2  <°»  GA_Gmin^  (2.66) 

where  G  =  4g  and  p  =  .3  is  chosen  to  match  the  Cf  (t)  =  a  , (t)  data.  The 
min  g  er  yl 

rationale  for  this  assumption  is  as  before:  under  high  G-stress  the  subject  tends 
to  tolerate  larger  tracking  errors.  These,  in  turn,  result  in  a  lower  relative 
weighting  on  the  error  in  the  cost  functional. 


Observation  Noise:  This  parameter  was  applied  precisely  in  the  same  manner  as  in 
the  longitudinal  case: 


Vyi(t)  =  "pyi 


fi*Ni 


i-1,2,3,4 


(2.67) 


The  nominal  p  values  are  -21.5  dB  as  before,  and  f ^  *  1/4.  Notice,  however,  that 
the  yi  term  is  absent  from  Equation  (2.66)  as  there  is  no  target  input  excitition, 
and,  strictly  speaking,  y^  *  0. 


Motor  Noise:  In  the  Optimal  Control  Modeling  process,  the  control  6g(t)  is  generated 
as 

T„  .  6  ’  ( t )  +  6  (t)  =  L ' (t)  x(t)  +  v  (t).  (2.68) 

N, equiv  a  a.  -  v  —  u 

The  gains  L(t)  are  computed  via 

L(t)  *  L*/K(t)  (2.69) 

where  L*  are  the  optimal  gains  corresponding  to  K(t)  -  1.  The  quantity  v^(t)  is  the 
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human's  motor  noise.  The  motor  noise  consists  of  two  parts:  a  purely  additive  com¬ 


ponent  v°(t)  and  a  multiplicative  component  that  scales  with  the  control 
u 

Thus, 


*'(t). 

a 


v  (t)  =  p  ( t )  v  °(t)  +  f,  ( t )  1 6  (t)| 
u  a  u  a 

where  :(t)  and  v  °(t)  are  assumed  to  be  white  and 
u 

modifies  the  additive  motor  ncise  covariance  (see 
Equation  (2.68)  through  by  K(t) ,  noting  Equation 
we  obtain 


(2.70) 

p  (t)  is  a  known  time  function  that 
a 

subsequent  discussion).  Multiplying 
(2.69)  and  defining  Sa(t)K(t)=u(t) , 


1  .  u(t)  +u(t) 

N, equtv 


x 

— o 


* 

-L 


=  A  x  +  b  u(t) 
o—o  — o 

'x(t)  +  Pa(t)K(t)v°(t)  +  «t)|u(t) 


(2.71) 

(2.72) 


Thus,  tiie  only  effect  of  the  gain  K(t)  is  to  increase  the  effective  additive  motor 

2 

noise  component,  the  variance  of  which  increases  with  K  (t).  The  motor  noise 
covariance  becomes  therefore 


V  (t)=V°(t)+iTp  ( t)  cov  [u(  t)  ]=p^  (t)K^  (t  )V°-Htp  (t)cov[u(t)]  (2.73) 

u  u  u  a  u  u 

where  p  (t)  is  the  motor  noise/signal  ratio,  and  irp^t)  is  identified  as  the  covariance 
of  f,(t).  The  multiplicative  part  of  the  motor-noise  scales  with  the  effective  control 
u(t),  which  is  convenient  for  modeling.  The  net  result  is  as  expected:  A  constant 
motor  noise  input  will  have  a  greater  effect  on  lateral  tracking  error  during  periods 
of  high  pitch  rate  since  the  system  sensitivity,  K(t),  increases  as  q  =  y  increases. 

In  addition,  one  has  to  consider  the  crossfeed  effect  in  the  manipulator.  It 

^Notice  that  an  additive  component  must  be  assumed  here,  as  there  would  be 
no  other  excitation  to  the  system. 
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appears  that  during  the  "active"  target  maneuver  intervals  (i.e.,  when  G  =  z(t)^O), 

the  additive  motor  noise  should  increase.  During  these  tracking  periods,  the  pilot 
*  c  2 

commands  C,^  '  .75  g/sec  (see  Figure  27,  t>.[10.,  17.14]  sec  and  t  [27.)4,  31.14]  sec). 

Because  of  this  longitudinal  control-rate,  some  additive  noise  is  injected  into  the 

lateral  axis,  introducing  random  lateral  errors. 

The  above  discussion  is  relative  to  the  interdependencies  between  the  system 

dynamics  and  the  motor  noise.  The  effect  of  C-stress  on  this  motor  noise  remained 

to  be  determined.  In  the  course  of  the  modeling  work,  it  became  apparent  that  model 

predictions  were  quite  sensitive  to  variations  in  p  and  p  .  Increases  in  these  para- 

a  u 

meters  greatly  increased  the  magnitude  of  the  model-predicted  aej.(t).  Large  increases 
in  o  (t)  are  observed  between  the  static  and  the  dynamic-G  conditions  (Figures  24, 

26).  We  postulate  therefore  that  G  and/or  G  stress  increases  the  motor-noise.  In 
a  previous  effort  to  model  pilot  performance  under  vibration  stress,  Levison  [44] 
also  needed  to  increase  the  motor  noise  to  match  the  data.  Although  vibration  is  a 
form  of  physical  stress  different  from  sustained  acceleration,  it  is  possible  to 
extrapolate  from  Levison 's  results  to  our  study.  The  following  approach  was  taken. 


1.  The  nominal  value  of  =  -18,0  dB.  In  the  dynamic  conditions  we  select 
P^  =  -16.5  in  accordance  with  the  longitudinal  tracking  task  model. 

2.  In  the  static  condition  the  nominal  value  of  V°  =  1.  (rad/sec)^  was  chosen 

on  the  basis  of  the  manipulator  dead-zone  characteristics.  In  the  dynamic 

conditions  V°  =  2.  The  value  of  p  is  nominally  selected  as  1. 
u  a 

3.  Because  of  the  crossfeed  phenomenon  we  increase  the  motor  noises  at  the  times 
for  which  z(t)jt0.  This  is  illustrated  in  Figure  47.  Basically,  remains 
unchanged  in  the  static  condition  but  it  increases  to  -13  dB  under  G-stress. 
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Hq.  47  MOTOR  NOISES  IN  THE  LATERAL  TRACKING  TASK 


In  the  static  condition  wo  set  p  =  3  and  in  the  dynamic  p  =  10.  These  increases 

cl  3 

reflect  nine  -  and  1()0  -  times  higher  levels  of  V°(t)( respectively  (see  Equations  2.72 
and  rhis  indicates,  perhaps,  that  sustained  acceleration  stress  enhances  the 

manipulator  cross feed  effect. 

The  approach  taken  here  modifies  somewhat  the  prelimiary  model  of  (45). 


V 


.V. '.•<?*  ”i<  ‘  n*  T f-nte  •  The  model  was  found  to  be  <|uite  insensitive  to  variations 

in  the  neuromuscular  time  constant,  i  .  The  effective  time  constant  used  in  this 
studv  was 

t  =  T  +  t.  =  .3  sec. 

N.equiv.  N  t 


i\4'jm_znoc  Index:  The  optimal  control  policy  is  that  which  minimizes  the  cost 
f unct ional 


T  ( u) 


=  £im  q:  F. 

X*wo  T 


lqeNl(al)er(t) 


q^2(t) 


KU2(t) | 


dt 


(  2  .  7  A ) 


This  is  in  accordance  with  the  assumption  that  the  subject  wc  ighs  the  perceived,  rather 
than  the  displayed,  error.  The  weightings  that  were  chosen  to  obtain  the  optimal 
control  gains  were  q^  =  .01,  =  .1  in  the  static  condition,  and  qg  =  .01,  q^  =  .25 

in  the  dynamic  one.  Notice  the  much  higher  j(t)  weighting  which  indicates  the  higher 
relative  importance  of  the  roll  angle  in  the  control.  Also,  it  is  evident  that  under 
G-stress  the  relative  weighting  on  er(t)  is  lower  than  in  the  static  condition  (.01:. 25 
vs  .01:.!).  The  net  effect  is  similar  to  that  of  increasing  the  indifference  thres¬ 
hold:  both  decrease  the  control  gain  on  the  error,  thus  generating  higher  o  (t). 


2 .4.)  Modeling  Results  of  the  Lateral  Tracking  Task 


An  excellent  indication  of  performance  degradation  under  G-stress  is  the  higher 

magnitude  of  a  (t).  We  observe  that  o  (t)  in  the  Dynamic- C  condition  is  about  60 
er  er 

millradiands  during  the  peak  G  period,  whereas  in  the  Static  condition  it  is  about 
30  mrads.  The  model  reproduces  these  trends;  the  modeling  results  are  illustrated  in 


(■inures  4S  31.  (All  figures  relate  to  the  first  60-second  period  of  the  1-peak 
target  pro  file.) 

•  I  igwrt  4«:  SD  of  tracking  error,  Static  G. 

•  i  ii'ur,  49:  SD  of  tracking  error,  Dynamic-G. 

•  figures  jo  -  5  I  illustrate  the  model-vs-data  comparisons  of  o  (t)  under 

er 

the  two  conditions. 

2.9.4  Discussion 

It  is  quite  evident  that  the  model  reproduces  the  experimental  data  very  well. 
It  should  be  emphasized  that  no  modifications  to  the  OCM  nominal  parameters  were 
necessary  in  the  static  condition.  This  reaffirms  the  validity  of  the  OCM  when 
nominal  conditions  are  considered.  Similarly  to  the  longitudinal  task,  the  increase 
in  the  model-predicted  o  (t)  in  the  dynamic  condition  is  attributed  mainly  to  the 
increase  in  the  human  motor-noise,  the  lower  error  weighting  in  the  cost  functional, 
and  the  increase  in  the  tracking  error  threshold  parameter. 
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Fig.  49  STANDARD  DEVIATION  HEADING  ERROR,  1  PEAK,  DYNAMIC  G 
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III.  FREQUENCY  DOMAIN  MEASURES  OF  HUMAN  PERFORMAMCK  I 'N I )HR  J ' -  ST K KS  ; .  _ ' 

In  Chapter  II,  the  development  of  a  "Ti  me-boma  in"  0— node  l  is  discussed.  Although  this 
performance  model  issuccessfulinthesen.se  that  it  was  validated  bv  the  experimental 
data,  the  information  it  furnishes  is  somewhat  limited.  It  is  intrinsically  incapable 
of  providing  any  frequency  measures,  such  as  power  spectra,  transfer  functions,  rem¬ 
nants.  In  fact,  this  model  is  only  effective  in  predicting  pilot  performance  in 

the  presence  of  transients  in  the  target  profiLe;  steady-state  analysis  is  excluded. 

It  was  imperative,  therefore,  to  build  a  stationary  performance  model  that  will 
complement  the  nonstationary  one.  In  this  chapter,  a  steady-stale  model  or  pilot 
performance  under  constant  sustained-G  stress  is  presented.  The  experimental  program 
leading  to  this  G-model  development  is  still  in  progress  at  AMRL.  Results  of  the 
initial  l)ES  runs  have  become  available  and  they  comprise  the  data  bast  for  the 
present  model  building  effort.  Three  experimental  conditions  are  employed,  and 
are  designated  as  (1)  Static  (fixed-base),  (2/  Semi -dynamic  .5  Cl  open-loop  stress1  a:id 
(3)  Post-runs  (static  runs  immediately  following  the  semi-dynamic,  condition).  The 
control  task  is  compensatory  tracking  in  the  lateral  plane.  The  target  input  excita¬ 
tion  is  a  psuedo-random  signal  (sum  of  sine  waves),  providing  a  rich  set  of  measures 
amenable  to  FFT  analysis.  This  is  discussed  in  detail  in  section  3.1.  The  perfor¬ 
mance  model  is  based  upon  the  Optimal  Control  Model  (OCM).  Section  3.2  describes  the 
modeling  effort. 

The  experimental  data  exhibit  significant  tracking  performance  decrements  under 
G-stress.  The  pilot's  control  remnant  increases  by  as  much  as  5dbai  some  frequencies, 
and  his  tracking  RMS  error  increases  by  some  15%.  In  addition,  the  pilot's  describing 
function  magnitude  decreases  by  nearly  1.5 db,  particularly  in  the  post-runs  condition. 

^This  Chapter  is  based  upon  a  1979  CDC  paper  1 48). 
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The  model  reproduces  these  trends  quite  well,  by  properly  adjusting  its  parameters. 
A  discussion  of  the  empirical  and  modeling  results  follows  in  Section  3.3. 


3 . _1  E xperimental  Design 

The  constant-G^  experiments  that  provide  the  data  base  for  the  modeling  effort 

reported  herein  were  not  designed  as  experiments  per  se.  In  fact,  they  are 

embodied  in  a  larger  experimental  program  designed  to  investigate  the  combined  G  -G 

z  y 

effects  on  pilots  in  the  new  Vectored  Force  Fighter  (VFF)  aircraft.  The  AMRL  DES 
centrifuge  was  the  motion  simulator  used  in  these  experiments.  The  subjects  were 
seated  in  the  cab,  which  represents  an  aircraft  cockpit,  and  was  equipped  with  a 
sidearm  control  stick.  The  visual  tracking  loop  included  a  compensatory  tracking  task  in 
the  lateral  plane  as  shown  in  Figure  52.  The  display  included  a  stationary  gunsight 
reticle  and  a  target  aircraft,  represented  by  a  dot. 


CONTROLLED  VEHICLE 

NOISE  SHAPING 

DYNAMICS 

FILTER 

1.7 

S  +  1.7 


+  s  +  .5 


CONTROL 

u(t) 


DISPLAY 


MANIPULATOR 


n(t)  -  target  excitation  signal 
8  (t)  -  vehicle  pointing  angle 


wQ(t)  -  white  driving  noise 


e(t)  -  tracking  error  =  n(t)  -  6D(t) 


52  TRACKING  TASK  CONFIGURATION 


and  solve  the  covariance  equation 


X  a'  +  A  X  +  E  WE'  =  0,  X  =  E{x  x'} 
s  s  s  s  — s  — s  s  — s— s 


(3.4) 


we  obtain  W  =  .36,  which  corresponds  to  the  specified  O^. 

The  random  signal  n(t)  was  simulated  by  combining  11  sine  waves  into  psuedo- 

random  forcing  function, 

11 


n(t) 


^  ^  A^  cos(w^t  +  0^) 


(3.5) 


for  each  run,  the  phases,  0^,  of  the  sinusoids  were  randomized  in  order  to  eliminate 
cognitive  response,  and  the  frequencies,  were  chosen  to  span  the  range  of  the 
human  response  (.15  :  15.57  rad/sec).  The  amplitudes,  A^ ,  were  selected  to  approxi¬ 
mate  |>f^(uj).  A  summary  of  the  frequency  components  comprising  the  pseudo-random 
signal  is  given  in  Table  1.  A  detailed  description  of  how  to  generate  the  approxi¬ 
mated  amplitudes  can  be  found  elsewhere  [46). 


No. 

Frequency  (rad/sec) 

Amplitude 

1 

0.15 

0.2376 

2 

0.38 

0.2404 

3 

0.69 

0.2023 

4 

1.00 

0.1330 

5 

1.46 

0.0854 

6 

2.22 

0.0482 

7 

3.  30 

0.0259 

8 

4.76 

0.0147 

9 

7.29 

0.0081 

10 

10.51 

0.0044 

11 

15.57 

_ 

0.0025 

Table  1  SUMMARY  OF  FREQUENCY  COMPONENTS 
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The  VFF  is  designed  in  such  a  way  that  the  pilot's  lateral  commands  generate  a 

transient  direct  side  force,  i.e.,  a  0^  acceleration.  In  order  to  simulate  the  desired 

t>y  levels,  it  was  necessary  to  centrifuge  the  subjects  to  an  (almost)  constant  G , 

and  to  rotate  the  DES  cab  (in  conjunction  with  the  pilot  tracking  commands)  to  produce 

an  off-normal  G  component.  The  present  study  deals  only  with  the  constant  G  runs 

z 

(5 G level)  and  no  G  acceleration.  The  tracking  task  was  designed  to  be  used  both  in 
the  present  constant  G^  study  as  well  as  in  future  (closed-loop)  G  applications. 

This  design  constrained  the  driving  signal  and  the  plant  dynamics  to  parameters 
comensurate  with  the  DES  cab's  bandwidth.  The  design  specifics  are  listed  below. 

Target  Input  Exaltation  Signal:  The  target  driving  signal  was  chosen  to  be  a  second 

order  colored  noise  n(t)  with  the  power  spectral  density 


VW)  = 


,  25W 


(3.1) 


a)  +  .25 


and  with  the  prescribed  RMS  value  o  f  a  =jEin2}  =  .3.  Such  noise  is  normally  generated 

as  the  output  of  a  linear  time  invariant  system,  T(s),  driven  by  white  noise  w^(t) 

2 

with  E{wo(t)}  =  W.  For  the  given  $^^(w),  this  system  is 


T(s)  = 


.5 


(3.2) 


s  +  s  +  .  5 

which  is  nothing  but  a  Butterworth  filter  with  a  cutoff  frequency  of  .707  rad/sec. 
If  we  represent  T(s)  in  the  state  space. 


V-s  +  Wt)= 


0 

1 

' 

.  X  + 

'  o  ' 

-.5 

-1 

— s 

.5 

wo(t) 


(3.3) 


H(t)  =  H  x 
-  s— s 


[1  0] 


X 

— s 
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.'antra  led  Vehicle  Dynamos:  The  chosen  plant  had  a  second-order,  overdamped  trans¬ 
fer  function,  which  was  dictated  by  the  DES  cab  dynamics,  viz.. 


C(s) 


3.4 

(s+1) (s+1. 7) 


(3.6) 


Run  Lenjth:  In  order  to  comply  with  medical  safety  requirements,  the  run  duration 
had  to  be  limited.  A  run  length  of  81.92  seconds  was  chosen;  with  a  sampling  period 
of  40  msec,  each  run  resulted  in  2048  datum  points,  facilitating  the  FFT  analysis. 

Inde:  ancient  Variables :  There  were  three  experimental  conditions  (excluding  the 

Dynamic  G  case), 
y 

1.  "Static"  -  The  centrifuge  cab  served  as  a  fixed-base  simulator  (no  Gz~stress) 
The  subject's  stick  inputs  effected  a  displayed  tracking  error  with  no 
actual  acceleration. 

2.  "Semi -dynamic"  -  The  centrifuge  cab  served  as  a  moving-base  simulator.  The 
pilot  nominal  level  was  5  G.  As  in  1,  the  subject's  control  corrected  the 
displayed  tracking  error,  but  no  off-normal  accelerations  were  present. 

3.  "Post-runs"  -  Same  as  1,  but  these  runs  immediately  followed  condition  2. 

Each  session  included  two  static  runs,  followed  by  two  semi-dynamic  runs, 
immediately  followed  by  two  post-runs.  By  this  arrangement  it  was  hoped  to  uncover 
some  of  the  sustained-,  as  well  as  the  accumulative-G  effects  on  tracking  performance 


Re: lie  it  ions:  Five  subjects  participated  in  a  total  of  nine  sessions.  With  two  runs 
per  condition  per  session,  18  runs  per  condition  resulted. 
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Data  Analysis:  Numerous  quantities,  ranging  from  tracking  error  to  heart  rate  time 
histories  were  recorded,  but  only  a  subset  of  these  quantities  was  required  for  the 
performance  model.  The  measurements  that  are  of  consequence  to  the  modeling  effort 
were  the  tracking  error,  e(t),  and  the  pilot  stick  command,  u(t). 

When  analysing  the  experimental  data,  the  following  steps  were  taken: 

1.  The  raw  data  of  each  of  the  indi”'dual  runs  were  transformed  to  frequency 

domain  via  an  FFT  program,  which  produced  the  following  quantities  (evaluated 
at  the  input  frequencies): 

(a)  Ensembles  of  the  individual  power  spectra  of  the  tracking  error,  e(t), 
and  of  the  pilot  control,  u(  t)  ($^^(01)  and  ^(u),  respectively). 

(b)  The  spectra  1>ee((jj)  and  $^(01)  are  given  in  two  parts,  the  input  - 
correlated  and  the  input  -  uncorrelated: 


<t  (u>)  =  <J>C  (w)  +  <t>uc  (to) 
ee  ee  ee 


<t>  (w) 

uu 


=  <tC  (00)  +  <!>UC(u>) 

uu  uu 


(3.7) 


(c)  The  cross-power-spectral  densities  of  e(t)  and  of  u(t)  with  the  input 
signal  n(t)  (*  (ju>)  and  ^ ( j uj) ,  respectively)  were  computed,  and  an 

ensemble  of  the  equivalent  pilot  describing  functions,  H(jw),  (magnitude 
and  phase)  was  obtained 


H(jw)  = 


Vjai) 

Vjw) 


(3.8) 


(d)  The  pilot  control  remnant,  R(w),  was  next  obtained  from  the  ratio  of 


4>u^(w)  to  <T  (u>), 
uu  uu 


R(U>)  = 


.uc ,  , 
$  (w) 

uu 

.c 

<P  0) 


(3.9) 
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This  quantity  is  crucial  in  any  performance-modeling  efforts.  As 
performance  degrades  under  G-stress,  it  is  reasonable  to  expect  an 
increase  in  the  remnant's  magnitude,  indicating  perhaps  involuntary 
motor  activity.  This,  indeed,  is  the  case  in  the  present  work,  as 
seen  in  the  subsequent  sections. 

(e)  In  addition,  the  performance  "scores,"  i.e.,  the  RMS  values  of  e(t) 

and  of  u(t)  (or  alternately,  the  total  power  in  1>ee  (<*>),  and  4>uu(w)  were 
computed.  These  quantities  prove  to  be  indicative  of  pilot-tracking 
performance,  as  they  increase  under  G-stress,  as  might  be  expected. 

2.  The  data  were  ensemble  averaged,  utilizing  techniques  similar  to  those  of 

Section  2.2.  The  final  product  was  first-  and  second-order  statistics  of: 

a.  Pilot  describing  function,  H(jui)  (magnitude  and  phase) , 

b.  Pilot  remnant,  R(u>), 

c.  RMS  scores. 

These  results  are  illustrated  in  figures  53-61  and  in  table  2. 

Discussion  of  the  experimental,  as  well  as  the  modeling  results  follows  in 

Section  3.3. 

3.2  Modeling  Approach 

The  base-line  model  adopted  here  was  the  "Steady-State"  OCM,  as  explained  in 
Section  1.4. 2. 2.  The  system  being  controlled  is  described  by  the  state  space  equation 


x  (t)  =  A  x  (t)  +  b  u(t)  +  E  w  (t) 

— o  o—o  — o  -  -o  o 


(3.10) 


V 


’ 

0 

1 

0 

0 

0 

0 

-1 

0 

0 

x  (t)  + 

0 

u(t)  + 

.5 

0 

0 

-1.7 

1  .  7 

-~o 

0 

0 

i  0 

0 

0 

-1 

2 

0 

1 

- 

- 

whori-  tho  tirst  two  states  are  associated  with  the  input  noise  shaping,  u(t)  is  the 

pilot’s  control  input  and  w  (t)  is  a  zero-mean  white  daussian  noise,  with 

o 

K  i w  ^ ( t )  ( ■ ) )  =  W6 (t-o) ,  W  =  . 36  (3.11) 

Tlu*  displayed  variables  consist  of  the  paired  tracking  error/error  rate,  viz.. 


1 

r  l 

e 

10-1  0 

=  C  x  (t)  = 

* 

0—0 

0  1  1.7  -1.7 

L  - 

- 

(3.12) 


The  OCM  assumes  that  the  well-trained  and  motivated  human  operator  adopts  an 
optimal  control  strategy,  subject  to  his  inherent  psycho-physiological  limitations. 
The  mathematical  interpretation  of  this  statement  is  that  the  pilot,  in  the  case  at 
hand,  will  minimize  the  following  cost  functional: 


J(u)  =  E  |  £im 
I  T-*°° 


if 


(e2(t) 


gu2(t] 


J 


(3.13) 


where  g  is  the  relative  weighting  on  the  control  rate.  Obtaining  the  optimal  control 
law  is  accomplished  by  solving  the  steady-state  Riccati  equation  pertaining  to  (3.13). 
The  well  known  solution  to  this  problem  is 


u+u  =  -L 


-opt 


x  (t) 
-o 


V  (t) 
u 


(3.14) 


i 
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whore  r  is  the  human  operator  (HO)  neuromuscular  time  constant,  v  (t)  is  an  additive 
N  u 

motor  noise,  L  is  the  optimal  control  gains  vector  and  x  (t)  is  the  estimated 
—  opt  x> 

state.  There  is  a  one-to-one  correspondence  between  q  and  r  usually  the  latter  is 

r  N 

specified  and  q  is  adjusted  accordingly. 

Human  Limitations :  The  neuromuscular  time  constant,  t  .  is  just  one  of  the  OCM  para- 

N 

meters  that  reflect  the  human  inherent  constraints.  It  represents  the  human's  neuro¬ 
motor  bandwidth  limitation  in  generating  the  control  input. 

Another  human  constraint  stems  from  the  fact  that  the  human  estimates  the  system 
state  on  the  basis  of  a  delayed  and  noisy  perceived  information, 

ypi ( t)  =  y.(t-x)  +  v  (t— i)  i=l,2  (3.15) 


where 


t  =  human's  lumped  time  delay, 

v  ,(t)  =  observation  noise  of  the  i-th  indicator, 
vi 

The  observation  (or  sensor)  noise  v^^(t)  is  a  zero-mean,  white  Gaussian  noise  with 
covariance 


E  {v  .(t)  v  .(a)}  =  V  ,<5(t -O) 
yi  yi  yi 

2 

V  .  =  TTp  .  E{y. (t) }, 
yx  yi  1 

and  p  is  the  noise/signal  ratio  of  the  i-th  observed  variable. 

The  driving  motor  noise  of  Equation  (3.14)  is  also  zero-mean. 
Its  covariance  is 


(3.16) 


white  and  Gaussian. 
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(3.17) 


E{v  (t)v  (o)}  =  V  6(t-CJ) 
u  u  u 

V  =  up  E{u2(t) }  , 
u  u  ’ 


where  p  is  the  motor  nosie/signal  ratio, 
u 

In  the  modeling  work,  r  ,  T,  p  and  p  are  of  ma]or  interest,  since  these  are 

N  y  i  u 

the  only  parameters  that  need  to  be  tuned  in  order  to  match  the  experimental  data. 

A  comprehensive  OCM  parameter  identification  program  that  was  developed  at  the 
University  of  Connecticut  [47]  facilitated  this  task.  The  application  of  this 
identification  procedure  to  the  case  at  hand  can  be  found  in  the  Appendix.  The 
experimental  and  the  modeling  results  are  discussed  in  the  following  section. 


3.3  Results  and  Discussion 

Kieures  53  and  61  present  the  model-data  comparisons  for  all  three  condi¬ 
tions.  The  data  are  represented  by  the  ensemble  means  and  the  standard  deviations 
at  the  input  frequencies,  and  the  model  predictions  —  by  the  continuous  lines.  In 
addition,  the  predicted  RMS  control  and  tracking  error  scores  are  compared  with  the 
experimental  ones  in  Table  2.  The  following  observations  can  be  made,  regarding 
the  empirical  data: 


Describing  Function  Magnitude:  The  highest  magnitude  among  the  three  conditions  is 
observed  in  the  static  case.  It  is  0.5  -  1.0  dB  larger  than  the  magnitude  in  the 
semi-dynamic  condition,  and  1.0  -  1.5  dB  larger  than  in  the  post-run  condition. 
These  differences  are  present  at  almost  all  input  frequencies. 


Describing  Function  Phase:  Here,  no  significant  differences  are  apparent  among  the 
three  conditions. 


102 


V 


Control  Remnant:  The  most  pronounced  difference  exists  here.  Whereas  there  is 
hardly  any  difference  in  the  remnant  between  the  static  and  the  post-run  conditions, 
the  semi-dynamic  case  (G^  -  stress.  Figure  3.7)  exhibits  an  increase  by  as  much  as 
)  to  5  d8  throughout  the  entire  frequency  range. 

Tmiakinj  Scores:  Table  2  lists  the  empirical  RMS  values  (+  one  standard  deviation) 
of  e(t)  and  of  u(t)  for  the  three  conditions.  The  following  observations  are  made: 

1.  The  lowest  RMS  error  value  is  in  the  static  condition,  and  the  highest--in 

the  G  -stress  condition, 
z 

2.  The  lowest  control  levels  were  applied  during  the  post-runs,  and  the  highest — 

during  the  G  runs, 
z 

"  i  .’ures  51  through  6  1  show  the  modeling  results,  which  are  the  major  achievement 
of  this  effort.  The  experimental  data  validate  the  model  in  all  three  cases,  through¬ 
out  the  entire  frequency  range.  Some  discrepancies  exist  at  low  frequencies,  between 
the  phases  in  particular:  the  model  is  unable  to  reproduce  the  phase  droops  below 
0.7  rad/ sec.  This  has  traditionally  been  the  case  in  all  past  modering  works. 

The  model  predicted  scores  match  the  experimental  RMS  values  as  well.  The  pre¬ 
dicted  scores,  both  error  and  control,  deviate  no  more  than  half  a  standard  devia¬ 
tion  from  the  ensemble  score  means. 

Table  3  summarizes  the  OCM  parameters  that  were  identified  from  the  data.  It 

is  apparent  that  the  operator's  motor  mechanism  is  impaired  by  Gz-stress.  The 

highest  motor  noise/signal  ratio,  p  ,  was  identified  for  the  semi-dynamic  condition. 

The  increase  in  p  under  G-stress  is  bv  as  much  as  A. 5  dB.  This  result  is  in  excellent 
u 

agreement  with  past  psycho-motor  research  [13,  19],  and  may  be  attributed  to  a  phenomenon 
such  as  tremor.  It  is  quite  interesting,  however,  that  the  post-runs  effect  lower 
mot<-  r  noise. 


1 03 
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Moreover,  the  observation  noise/signal  coefficient,  p  ,  under  G-stress  increases 

from  -19.6  dB  to  a  level  as  high  as  -15.1  dB.  This  again  is  consistent  with  past 

research  concerning  G-effects  on  vision.  The  increase  mav  be  accounted  for  by 

decrements  in  visual  acuity  under  G-stress,  blur,  etc.,  as  reported  by  (aether  [h; 

and  manv  others.  Alternately,  this  increase  in  i  can  be  ascribed  to  a  decrease  in 
■  '  e 

the  attentional  allocation  or  to  an  increase  in  the  indifference  threshold  on  the 
observed  error  due  to  G-stress,  or  both.  This  is  readily  evident  i f  we  replace 
Equation  (3.1b)  with 


Tip  E* e  (t)  . 
e _ _ 

f  -N 
e  e 


(3.18) 


where  f  <  1  is  the  fractional  attention  allocated  to  e(t).  The  term 


N  =  N  (a  ,  E{e  (t ) }  ]  <  1 
e  e  e  — 


(3.19) 


is  a  statistically  linearized  function,  arising  from  a  thresholded  error  observation 
with  threshold  level  ag  (see  (43]).  It  has  been  shown  [44,  45]  that  the  error  in¬ 
difference  threshold  is  indeed  affected  by  G-stress. 

The  trends  in  p*  are  not  as  clear  as  in  p  ,  but  the  model  proved  quite  insensitive 
e  e 

to  the  former.  Therefore,  may  be  assigned  a  prescribed  nominal  value. 

Another  parameter  which  is  affected  bv  G-stress  is  the  lumped  time  delay,  t. 

There  is  some  increase  in  r  in  the  semi-dynamic  condition  but  it  decreases  slightly 

in  the  post-run  case.  The  increase  in  human  reaction  time  under  G-stress  is  again  in 

accordance  with  previous  research  (e.g.  [13]).  The  neuromuscular  time  constant,  t^, 

increases  both  for  G  and  post-G  runs;  this  validates  the  hypothesis  of  accumulative 

G  -effects.  This  decrease  in  the  human  control  bandwidth,  combined  with  a  significant 
z 
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decrease  both  In  the  describing  function  gain  (Figure  59)  and  in  control  activity, 
(Table  2  -  post  run  condition),  suggest  a  hystersis  effect  of  fatigue  induced  by 


an  accumulated  G  -stress. 

z 

I'ABLF  2  -  RMS  SCOKKS 


SCORES 

^(RMS) 

CONDITION 

TRACKING  ERROR,  e 

CONTROL 

u 

EXPERIMENTAL 

PREDICTED 

EXPERIMENTAL 

PREDICTED 

STAT 1 C 

.093  +  .014 

.090 

.37  +  .10 

.32 

G  -STRESS 
z 

.107  +  .015 

.  102 

.42  +  .18 

.  33 

j  POST-RUNS 

.095  +  .021 

_ 

.088 

.33  +  .09 

.29 

Table  3  -  MODEL  PARAMETERS 


^CONDITION 

parameters'^ 

STATIC 

G  -  STRESS 
z 

POST-RUNS 

P  y 

U 

dB 

-17.5 

-15.3 

-20.1 

Ue  * 

dB 

-19.6 

-15.1 

-21.3 

• 

dB 

-29.5 

LO 

O 

00 

-24.5 

T  , 

SEC 

.257 

.266 

.2  36 

tn  ’ 

SEC 

.159 

.176 

.165 

Iu5 


1 


DEGREES 


>8 


DEGREE 


1 V  CONCLUSIONS 

It  has  become  a  well-established  fact  that  linear  acceleration  stress  effects 
substantial  degradation  in  pilot  performance.  The  results  of  the  compensatory 
tracking/sustained  acceleration  research  program,  jointly  undertaken  by  Air  Force 
Aerospace  Medical  Research  Laboratory  and  the  University  of  Connecticut -CYBERLAB, 
substantiates  this  contention.  Statistical  analyses  performed  on  the  experimental 
data  prove  that  significant  differences  exist  between  versus  no-G^-stress  condi¬ 
tions.  The  acceleration  stress  data  exhibit  trends  which  are  clearly  a  manifestation 
of  performance  impairment: 

1.  The  variability  in  the  tracking  errors  increases  significantly  under  Gz~stress. 
This  is  demonstrated  by  high  levels  and  enduring  periods  of  tracking-error- 
variance  time  histories,  and  by  larger  RMS  tracking  scores. 

2.  Pilot's  involuntary  motor  activity  increases  in  the  presence  of  sustained 
C^-stress.  The  data  show  large  increase  in  the  input-uncorrelated  control 
power,  and  higher  control  RMS  values. 

5.  it  is  evident  from  the  data  that  G  -stress  effect  slow  tracking  error 

z 

recovery  and  lower  magnitude  of  the  human  describing  function.  I  his  is, 
perhaps,  an  indication  of  a  lower  control  bandwidth  on  the  part  of  the  pilot 
under  positive  acceleration  stress. 

The  analytical  phase,  on  its  part,  has  resulted  in  a  normative  performance  model 
which  has  been  validated  by  the  experimental  data.  The  model  is  based  upon  the 
Optimal  Control  Model  (OCM)  and  it  has  been  augmented  by  a  G-submodel  to  include  the 
acceleration  effects.  The  model  parameters,  which  reflect  a  person's  inherent  limita¬ 
tions,  depart  from  their  nominal  values  to  account  for  the  performance  degradation 
under  strenuous  acceleration  conditions. 
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Virtually  all  the  model  parameters  need  to  be  perturbed  in  order  to  replicate  the 
data  trends.  It  is  shown  that  the  model  parameters  which  are  affected  are  the  observa¬ 
tion  noise,  neuiomotor  time  constant  t  ,  and  reaction  time-delay,  all  of  which  exhibit 

a  significant  increase  under  G  -stress.  In  addition,  some  accumulative  G-effects  are 

z 

evident  as  remains  higher  than  nominal  in  the  post-G  runs,  reflecting  a  lower 
control  bandwidth,  and  perhaps  fatigue. 

The  most  noteworthy  deviations  from  the  nominal,  however,  are  the  increases  in 
the  motor  noise  and  in  the  error  indifference  threshold  parameters.  The  first 
results  in  higher  control  remnant  levels  and  tracking  error  variances,  which  are 
observed  from  the  empirical  data.  The  latter  effects  a  slow  tracking  -  error  recovery 
in  the  closed-loop  G-stress  condition. 

At  present,  we  are  pursuing  a  G^-stress  research  program  which  is  an  extension  of 
this  study.  The  G^  program,  in  addition  to  lateral  accelerations,  involves  positive 
vertical  acceleration  stress  of  up  to  5  G.  As  the  experimental  data  of  this  program 
become  available,  the  performance  model  under  Gz-stress  may  be  refined  further. 


APPENDIX 

l DENT  l  FICATION  OF  OCM  PARAMKT ERE  UNDER  G- STRESS 

The  application  of  the  Optimal  Control  Model  (OCM)  to  a  manual-control  task 
usually  involves  matching  model  predictions  to  available  experimental  data.  Tradi¬ 
tionally,  this  has  been  done  bv  tediously  varying  the  model  parameters  until  a 
suitable  model-data  match  has  been  obtained.  In  the  past  year,  a  comprehensive 
Parameter  Identification  Program  for  the  OCM  (OCMID)  has  been  developed  [47).  This 
program  provides  an  automated  technology  for  the  model-building  process. 

The  OCMID  is  a  modular  computer  package  designed  to  identify  parameters  of  the 
OCM.  User-supplied  routines  determine  the  parameters  to  be  identified  and  the  form 
of  the  program  output.  The  package  is  designed  with  flexibility  to  allow  any  routine 
to  be  interchanged  with  user  -  supplied  routines,  so  the  package  can  be  tailored  to 
specific  needs.  In  conjunction  with  the  present  G-model  research,  the  OCMID  has  been 
extended  with  additional  routines,  incorporating  a  graphics  display  subprogram. 

Approach : 

The  objective  is  to  identify  the  pilot  describing  function  (magnitude  and  phase), 
the  remnant,  and  the  model  predicted  performance  scores.  The  OCM  parameters  that 
generate  these  predictions  are  to  be  identified  from  experimental  frequency  measures 
data  so  as  to  minimize  a  weighted  least-squares  objective  functional.  The  minimiza¬ 
tion  of  this  cost  functional  provides  a  scalar  metric,  which  in  turn,  reflects  how 
closely  the  currently  obtained  model  matches  the  measured  describing  function's 
magnitude  and  phase,  control  remnant  and  performance  RMS  scores.  We  choose,  therefore. 


M3 


the  cost  functional 


•  v,r,;;  ““ )  i 

sc  1  \  su  ' 


(A.l) 


N  .  i 

E'  I"  /  |IH  i'-^J  |“|Ht  i-iij  \  MU  |  .1 .  J  -  /H(  jo.  i  \  /  k  <:...)-R  (w.)\2" 

...l1"''  . ~’m  /  •  A . «pV .  )  V: . ) 


s  =  RMS  score  of  tracking  error 
=  RMS  score  of  control  input 
. )  =  Pilot  describing  function  i-th  input  fre<]uencv 
R  =  Pilot  remnant  (d  i-th  input  frequency 

'■T/.)  =  Ensemble  standard  deviation  of  the  pertinent  variable 
(•)  =  Represents  the  OCMID  current  estimates 
=  Relative  weightings  (at  present  =  1.) 

N^.  =  Number  of  frequencies  =  11 

The  choice  of  this  cost  functional  gives  a  weighted  least-squares  fit  to  the 
data.  Notice  that  each  squared  difference  in  .1  is  divided  by  the  appropriate  ensemble 
standard  deviation.  In  this  way,  datum  points  with  larger  a ^ ^  are  automatically 
given  less  weight. 

The  minimization  of  .1  is  carried  out  by  applying  Powell's  non-gradient  optimiza¬ 
tion  method  [49].  Basically,  Powell's  method  finds  pseudo-conjugate  gradient  search 
directions  by  perturbing  the  parameters  in  the  Euclidean  space.  A  bivariate  quadratic 
line  search  is  used  to  minimize  J  in  each  search  direction,  such  that  the  parameters 
"move"  a  step  along  each  orthogonal  direction.  This  means  that  in  one  Powell  itera¬ 
tion  the  parameters  "move"  along  the  average  direction  of  all  the  individual  orthogonal 
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di rect  ions . 
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After  the  program  has  converged,  a  single  quasi-Newton  step 
convergence  and  to  provide  approximate  confidence  intervals  for 
meters.  To  see  this  step  we  rewrite  J  as: 


is  taken  to  ensure 
the  identified  para- 


NT 

■I  =  y  Y.  f\  (6)  =  J 
i  =  1 


whe  re : 

NT  =  total  number  of  the  weigiited  squared  differences 

0  =  parameter  vector  to  be  identified 
2 

f^  =  i-th  squared  term  . 

The  quasi-Newton  step  is 


(A. 2) 


0k+1 


"  3f ' 

3f  ' 

_1  &  f 

30  ' 

30 

30- 

(A.  3) 


where  a  is  the  one-dimensional  line  search  minimizing  element.  In  the  case  at  hand. 


NT 

O' 


2  +  3-Nf  =  35 


(P. 


u. 


N’ 


T) 


and  the  f^(0)  terms  are  given  by  Equation  (A.l).  Note  that  the  term  inside  the  squared 
brackets  of  Equation  (A. 3)  is  the  approximated  Hessian.  The  approximate  confidence 
interval  for  the  identified  parameter  set  is  then  given  by  the  square-root  of  the 
diagonal  elements  of  the  inverse  Hessien. 
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